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SUMMARY 


An experimental and analytical investigation of the flutter of 
sweptback cantilever wings is reported. The experiments employed groups 
of wings swept back by rotating and by shearing. The angle of sweep 
ranged from 0° to 6 0° and Mach numbers extended to approximately 0 . 85 . 

A theoretical analysis of the air forces on an oscillating swept wing 
of high length— to— chord ratio is developed, and the approximations 
inherent in the assumptions are discussed. Comparison with experiment 
indicates that the analysis developed in the present paper is satis- 
factory for giving the main effects of sweep, at least for nearly 
uniform cantilever wings of high and moderate length— to— chord ratios. 

A separation of the effects of finite span and compressibility in their 
relation to sweep has not been made experimentally but some combined 
effects are given. A discussion of some of the experimental and theor- 
etical trends is given with the aid of several tables and figures. 


INTRODUCTION 


The present paper is an outgrowth of the trend toward the use of 
swept wings for high-speed flight and reports the results of an analysis 
and of an accompanying exploratory program of research in the Langley 
k. 5 — foot flutter research tunnel on swept cantilever wings. The 
material was assembled in a memorandum form with a similar title in 19 ^ 8 . 
The chief purposes of the present paper are to provide a more detailed 
exposition of the analysis and to make the main material more generally 
available . 

Mention of some previous experimental and analytical work on swept 
wings follows. A preliminary experimental investigation of the effect 
of sweep on flutter has been made (reference 1) with a single, simple 
rigid wing mounted flexibly at one end of a base which could be rotated 
to various desired sweep angles. This investigation was made at low 
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Mach numbers for two "bending— torsion frequency ratios and at several 
angles of sweepback. Another investigation (data unpublished) in which 
the density of the test medium was a variable was conducted by D. Benun 
on the same type of rigid, flexibly mounted wing at higher Mach numbers 
and at sweep angles of 0° and 45°. Other unpublished work on swept wings 
exists, but a search of the available information indicates a need for 
further systematic study. 

The experimental work reported herein dealt with models mounted as 
cantilevers at their roots. These cantilever models differed from the 
rigid, flexibly mounted wings, which had all bending and torsion flexi- 
bility concentrated at the root, and thus were subject to different root 
effects. In order to facilitate analysis the cantilever models were 
uniform and untapered. The intent of the experimental program was to 
establish trends and to indicate orders of magnitude of the various 
effects of sweep on flutter, rather than to isolate precisely the 
separate effects. 

The models were swept back in two basic manners — shearing and 
rotating. For the case in which the wings which were swept back by 
shearing the cross sections parallel to the air stream, the span and 
aspect ratio remained constant. For the other case, a series of 
rectangular— plan— form wings were mounted on a special base which could 
be rotated to provide any desired angle of sweepback. This rotatory 
base was also used to examine the critical speed of sweptforward wings. 

Tests were conducted also on special models that were of the 
"rotated" type (sections normal to the leading edge were the same at 
all sweep angles) with the difference that the bases were alined 
parallel to the air stream. Two series of such rotated models having 
different lengths were tested. 

Inasmuch as the location of the center of gravity, the mass— density 
ratio, and the Mach number have important effects on the flutter charac- 
teristics of unswept wings, these parameters were varied for swept wings. 
In order to investigate possible changes in flutter characteristics 
which might be due to different flow over the tips, various tip shapes 
were included in the experiments. 

In an analysis of flutter, vibrational characteristics are very 
significant; accordingly, vibration tests were made on each model. A 
special study of the change in frequency and mode shape with angle of 
sweep was made for a simple aluminum-alloy beam and is reported in 
appendix A. 

Theoretical analyses of the effect of sweep on flutter exist only 
in brief or preliminary forms. In England in 19^2 W. J. Duncan estimated 
by certain dimensional considerations the effect of sweep on the flutter 
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speed of certain specialized wing types. Among other British workers 
whose n a m e 3 are mentioned in connection with problems of flutter 
involving sweep are E. McEinnon Wood, A. E. Collar, and I. T. Minhinnick. 
An account of Minhinnick* s work was given by Broadbent in reference 2. 

In reference 3 a preliminary analysis for the flutter of swept wings in 
incompressible flow is developed on the basis of a "strip theory" (with 
the strips taken in the stream direction) and is applied to the experi- 
mental results of reference 1. Examination of the limiting case of 
infinite span discloses that the aerodynamic assumptions employed in 
reference 3 are not well-grounded. Eeference k adapts this "strip 
theory" to flexible wings and also presents an alternative "velocity 
component" treatment employing other aerodynamic assumptions which in 
their end result appear more akin to those employed in the analysis of 
the present paper. Wo definite choice is made in reference k between 
the two methods although the "strip theory" method is favored. 

In the present paper a theoretical analysis is developed anew and 
given a general presentation. Application of the analysis has been 
limited at this time chiefly to those calculations needed for comparison 
with experimental results. A wider examination of the effect of various 
parameters and of additional degrees of freedom on the flutter character- 
istics is desirable. 


SYMBOLS 


half— chord of wing measured perpendicular to elastic 
axis, feet 


r 


half— chord perpendicular to elastic axis at reference 
station, feet 


effective length of wing, measured along elastic axis, 
feet 


c 


wing chord measured perpendicular to elastic axis. Inches 


l 


length of wing measured along midchord line, inches 


A 


angle of sweep, positive for sweepback, degrees 


A 


g 


geometric aspect ratio 



coordinate perpendicular to elastic axis in plane of 
wing, feet 


coordinate along elastic axis, feet 
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z* 


Z 

*1 

I 

h 

9 


cr 


f h(y’), r h (n) 
fe(y*)^ F 0 (n) 


coordinate in direction perpendicular to x'y* plane , 
feet 

coordinate of wing surface in z* direction, feet 
nondimensional coordinate along elastic axis (y*/Z') 
coordinate in wind— stream direction 

bending deflection of elastic axis, positive downward 

torsional deflection of elastic axis, positive with 
leading edge up 

local 
local 

deflection function of wing in bending 
deflection function of wing in torsion 


bending slope of elastic axis 
rate of change of twist 



CD 

CDh 

f h, 


L hg 


time 

angular frequency of vibration, radians per second 
angular uncoupled bending frequency, radians per second * 

angular uncoupled torsional frequency about elastic axis, 
radians per second 

first bending natural frequency, cycles per second 
second bending natural frequency, cycles per second 
first torsion natural frequency, cycles per second 
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f 


a 


fe 


uncoupled first torsion frequency relative to elastic 


axis, cycles per second 




experimental flutter frequency, cycles per second 


f R reference flutter frequency, cycles per second 

f. flutter frequency determined "by analysis of present 

report, cycles per second 

v free— stream velocity, feet per second 

v e experimental flutter speed, feet per second 

v n component of air— stream velocity perpendicular to elastic 

axis, feet per second (v cos A) 

V e experimental flutter speed taken parallel to air stream, 

miles per hour 

V reference flutter speed, miles per hour 

rv 


reference flutter speed based on wing elastic axis, 
miles per hour (defined in appendix B) 

flutter speed determined by theory of present report, 
miles per hour 

theoretical divergence speed, miles per hour 

reduced frequency employing velocity component 
perpendicular to elastic axis 

cp phase difference between wing bending and wing torsion 

strains, degrees 

p density of testing medium at flutter, slugs per cubic 

foot 
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a 

M 

Mcr 


v ea 


ea 


a 


a + Xa 


m 


K 


I 


a 


r 


a 


El 

GJ 


& 


’h 


dynamic pressure at flutter, pounds per square foot 
Mach, number at flutter 
critical Mach number 


distance of center of gravity behind leading edge taken 
perpendicular to elastic axis, percent chord 

distance of elastic center of wing cross section behind 
leading edge taken perpendicular to elastic axis, 
percent chord 

distance of elastic axis of wing behind leading edge 
taken perpendicular to elastic axis, percent chord 


nondimens ional elastic— axis position 



nondime ns ional center— of— gravity position 



mass of wing per unit length, slugs per foot 


wing mass— density ratio at flutter 



mass moment of inertia of wing per unit length about 
elastic axis, slug— feet^ per foot 


nondimens ional radius of gyration of wing about elastic 



bending stiffness, pound— inches^ in tables, pound— feet^ 
in analysis 

torsional stiffness, pound— inches^ in tables, pound— 

o 

feet in analysis 


structural damping coefficient for bending vibration 
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structural damping coefficient for torsional vibration 


n 


a special bracket used to identify terms which are due 
solely to inclusion of the last term in equation (5a) 


Note: In order to preserve continuity and to facilitate comparison with 

previous work on the unswept wing, the subscript a rather than Q is 
retained with certain quantities to refer to the torsional degree of 
freedom. 


experimental investigation 

Apparatus 


Wind tunnel.— The tests were conducted in the Langley 4.5— foot 

flutter research tunnel which is of the closed-throat, single— return type 
employing either air or Freon— 12 as a testing medium at pressures 
varying from 4 inches of mercury to 30 inches of mercury. In Freon— 12, 
the speed of sound is 324 miles per hour and the density is 0.0106 slug 
per cubic foot at standard pressure and temperature. The maximum choking 
Mach number for these tests was approximately 0.92. The Reynolds number 

range was from 0.26 X 10^ to 2.6 x 10^ with most of the tests at Reynolds 
numbers of the order of 1.0 x 10^. 

Models .— In order to obtain structural parameters required for the 
flutter studies, different types of construction were used for the models. 
Some models were solid spruce, others were solid balsa, and many were 
combinations of balsa with various aluminumr-alloy inserts. Seven series 
of models were investigated, for which the cross sections and plan forms 
are shown in figure 1. 

Figure 1(a) shows the series of models which were swept back by 
shearing the cross sections parallel to the air stream. In order to 
obtain flutter with these low-aspect— ratio models, thin sections and 
relatively light and weak wood construction were employed. 

The series of rectangular-plan— form models shown in figure 1(b) were 
swept back by using a base mount that could be rotated to give the desired 
awee p angle. The same base mount was used for testing models at forward 
sweep angles. It is known that for forward sweep angles divergence is 
critical. In an attempt to separate the divergence and flutter speeds 
in the sweepforward tests, a L— spar cross-sectional construction was 
used to get the elastic axis relatively far forward (fig. 1(c)). 
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Two aeries of wings (figs. 1(d) and 1(e)) were swept back with the 
length—- to-chord ratios kept constant. In these series of models, the 
chord perpendicular to the leading edge was kept constant and the bases 
were alined parallel to the air stream. The wings of length— to-chord 
ratio 8.5 (fig. 1(d)) were cut down to get the wings of length— to— chord 
ratio 6.5 (fig. 1(e)). 

Another series of models obtained by using this same manner of 
sweep (fig. 1(f)) was used for investigating some effects of tip shape. 

Spanwise strips of lead were fastened to the models shown in fig- 
ure 1 (e) and a series of tests were conducted with these weighted models 
to determine the effect of center-of-gravity shift on the flutter speed 
of swept wings. The method of varying the center of gravity is shown in 
figure 1(g). In order to obtain data at zero sweep angle it was neces- 
sary, because of the proximity of flutter speed to wing-divergence speed, 
to use three different wings. These zero-sweep-angle wings, of 8 -inch 
chord and 48-inch length, had an internal weight system. 

The models were mounted from the top of the tunnel as cantilever 
beams with rigid bases (fig. 2). Wear the root of each model two sets 
of strain gages were fastened, one set for recording principally bending 
deformations and the other set for recording principally torsional 
deflections . 


Methods 

Determination of model parameters .— Pertinent geometric and struc- 
tural properties of the model are given in tables I to VII. Some 
parameters of interest are discussed in the following paragraphs. 

As an indication of the nearness to sonic— flow conditions, the 
critical Mach number is listed. This Mach number is determined by the 
Karman— Tsien method for a wing section normal to the leading edge at 
zero lift. 

The geometric aspect ratio of a wing is here defined as 

A = Semispan 9 = .cos a) 2 = 1 cog 2 A = A 

Plan— form area z c c 2 


The geometric aspect ratio Ag is used in place of the conventional 

aspect ratio A because the models were only semispan wings. For 
sheared swept wings, obtained from a given unswept wing, the geometric 
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aspect ratio is constant, whereas for the wings of constant length— to— 

chord ratio the geometric aspect ratio decreases with cos^\ as the 
angle -of sweep is increased. 

The weight, ce nt e r-of-gr av it y position, and polar moment of inertia 
of the models were determined by usual means. The models were stati- 
cally loaded at the tip to obtain the rigidities in torsion and bending, 
G J and El . 

A parameter occurring in the methods of analysis of this paper is 
the position of the elastic axis. A "section" elastic axis located 
a t x ga was obtained for wings from each series of models as follows: 

The wings were clamped at the root normal to the leading edge and at a 
chosen spanwise station were loaded at points lying in the chordwise 
direction. The point for which pure bending deflection occurred, with 
no twist in the plane normal to the leading edge, was determined. The 
same procedure was used for those wings which were clamped at the root, 
not normal, but at an angle to the leading edge. A different elastic 
axis designated the "wing" elastic axis and located at x * was thus 

6a 

determined. 

For these uniform, swept wings with fairly large length—' to-chord 
ratios, the "wing" elastic axis was reasonably straight and remained 
essentially parallel to the "section" elastic axis, although it was 
found to move farther behind the "section" elastic axis as the angle of 
sweep was increased. It is realized that in general for nonuniform 
wings, for example, wings with cut-outs or skewed clamping, a certain 
degree of cross stiffness exists and the concept of an elastic axis is 
an oversimplification. More general concepts such as those involving 
influence coefficients may be required. These more strict considera- 
tions, however, are not required here since the elastic— axis parameter 
is of fairly secondary importance. 

The wing mass— density ratio k is the ratio of the mass of a 
cylinder of testing medium, of a diameter equal to the chord of the 
wing, to the mass of the wing, both taken for unit length along the 
wing. The density of the testing medium when flutter occurred was used 
in the evaluation of k. 

Determination of the reference flutter speed .— It is convenient in 

presenting and comparing data of swept and unswept wings to employ a 
certain reference flutter speed. This reference flutter speed will 
serve to reduce variations in flutter characteristics which arise from 
changes in the various model parameters such as density and section 
properties not pertinent to the investigation. It thus aids n system- 
atizing the data and emphasizing the desired effects of sweep including 
effects of aspect ratio and Mach number. 
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This reference flutter speed may he obtained in the following 

way. Suppose the wing to be rotated about the intersection of the 
elastic axis with the root to a position of zero sweep. In this posi- 
tion the reference flutter speed is calculated by the method of refer- 
ence 5, which assumes an idealized, uniform, infinite wing mounted on 
springs in an incompressible medium. For nonuniform wings, a reference 
section taken at a representative spanwise position, or some integrated 
value, may be used. Since the wings used were uniform, any reference 
section will serve. The reference flutter speed may thus be considered 
a "section” reference flutter speed and parameters of a section normal 
to the leading edge are used in its calculation. This calculation 
also employs the uncoupled first bending and torsion frequencies of the 
wing (obtained from the measured frequencies) and the measured density 
of the testing medium at time of flutter. The calculation yields a 
corresponding reference flutter frequency which is useful in comparing 
the frequency data. For the sake of completeness a further discussion 
of the reference flutter speed is given in appendix B. 

Test procedure and records .— Since flutter is often a sudden and 
destructive phenomenon, coordinated test procedures were required. 

During each test, the tunnel speed was slowly raised until a speed was 
reached for which the amplitudes of oscillation of the model in bending 
and torsion increased rapidly while the frequencies in bending and tor— 
sion, as observed on the screen of the recording oscillograph, merged 
to the same value. At this instant, the tunnel conditions were recorded 
and an oscillograph record of the model deflections was taken. The tun- 
nel speed was immediately reduced in an effort to prevent destruction 
of the model. 

From the tunnel data, the experimental flutter speed V e , the den- 
sity of the testing medium p, and the Mach number M were determined. 
Wo blocking or wake corrections to the measured tunnel velocity were 
applied. 

From the oscillogram the experimental flutter frequency f e and 
the phase difference cp (or the phase difference ±l 80 °) between the 
bending and torsion deflections near the root were read. A reproduction 
of a typical oscillograph flutter record, which indicated the flutter to 
be a coupling of the wing bending and torsion degrees of freedom, is 
shown as figure 3* Since semispan wings mounted rigidly at the base 
were used, the flutter mode may be considered to correspond to the flut- 
ter of a complete wing having a very heavy fuselage at midspan, that is, 
to the symmetrical type. 

The natural frequencies of the models in bending and torsion at 
zero airspeed were recorded before and after each test in order to 
ascertain possible changes in structural characteristics. In most cases 
there were no appreciable changes in frequencies but there were some 
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reductions in stiffnesses for models which had been weakened by flut- 
tering violently. Analysis of the decay records of the natural fre- 
quencies indicated that the wing dam ping coefficients gh and go, 
(reference 5) were about 0.02 in the first bending mode and 0.03 in the 
torsion mode; 


ANALYTICAL INVESTIGATION 
General 


Assumptions.— An attempt is first made to point out the main 
assumptions which seem to be applicable for swept wings of moderate 
taper and of high or moderate length—' to-chord ratios. 

(a) The assumptions, such as small disturbances and potential flow, 
commonly employed in linearized treatment of unswept wings in an ideal 
incompressible fluid are made. 

(b) The structural behavior is such that over the main part of the 
wing the elastic axis may be considered straight. The wing is also con- 
sidered sufficiently stiff at the root so that it behaves as if it were 
clamped normal to the elastic axis. An effective length l* needed 
for integration reasons may be defined (for example, as in fig. 4 ). The 
angle of sweepback is measured in the plane of the wing from the direc- 
tion normal to the air stream to the elastic axis. All section parame- 
ters such as semichord, locations of elastic axis and center of gravity, 
radius of gyration, and so forth, are based on sections normal to the 
elastic axis. 

(c) The aerodynamic behavior is such that any section dy* of the 
wing normal to the elastic axis, taken in the direction of the compo- 
nent v cos A of the main— stream velocity, generates a velocity poten- 
tial associated with a uniform infinite swept wing having the same 
instantaneous distribution over the chord of velocity normal to the 
wing surface as does the actual section. 

Additional remarks on these assumptions are appropriate. With 
regard to assumption (a), in accordance with linearization of the prob- 
lem, the boundary conditions are stated and treated with respect to a 
reference surface, in this case a plane, containing the mean equilibrium 
position of the wing and the main-stream velocity. Furthermore, incom- 
pressible flow is assumed in order to avoid complexity of the analysis, 
although modifications due to iMach number effects can be added. Such 
modifications may be based, for example, for wings having large length— 
to— chord ratios, on existing theoretical calculations of aerodynamic 
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coefficients for subsonic or supersonic two-dimensional flow appropriate 
to the component v cos A. On the other hand the modifications may be 
partly empirical, especially for "transonic" conditions and for small 
length— to— chord ratios. The transonic conditions and the general aero- 
dynamic behavior of swept wings may depend, for large length— to— chord 
ratios, on the component v cos A, but the dependence may shift to the 
stream velocity v for small length—' to-chord ratios. 

With respect to assumption (b), results of analyses of and experi- 
ment on unswept wings having low ratios of bending frequency to torsion 
frequency show that small variations of position of the elastic axis 
are not important. The assumption of a straight elastic axis over the 
main part of a swept wing, similarly, is not critical for many cases. 
This assumption is made for convenience, however, and modifications for 
a curved elastic axis can be made when necessary, for example, for 
plate— like wings. Small differences in the angle of sweepback of the 
leading edge, quarter— chord line, elastic axis, and so forth, are neg- 
lected. The analysis could be further modified to take into account 
variation of the angle of sweepback along the length of the wing. 

Assumption (c) implies that associated with the action of the wing 
in pushing air downward there is a none ire ulatory potential— type flow 
similar to that around sections of an infinite flat— plate wing. 
Furthermore, as in the case of the unswept airfoil, a circulatory 
potential— type flow is generated in which for the swept airfoil the 
component v cos A is decisive in fixing the circulation. (This 
assumption differs from that made in the "strip theory" of references 3 
and k which employs the main— stream velocity together with sections of 
the wings parallel to the stream direction. ) Effects of the floating 
of the wake in the stream direction rather than in the direction 
of v cos A and induced effects of variation of the strength of the 
wake in the wing— length direction are neglected, as are three- 
dimensional tip effects. For large values of the reduced frequency k n 
a given segment of the wing might be influenced chiefly by the nearby 
wake and the correction would be small. On the other hand, for small 
values of k n a given segment might be influenced by a more widespread 
portion of the wake; corrections for this condition may possibly be 
based on knowledge of the static case (for example, slope of the lift 
curve). As the angle of sweep approaches 90 ° > obviously the mechanism 
for the generation of lift is different from the one postulated here; 
for example , a tip condition may replace the trailing— edge condition, 
and considerations of very small aspect ratio arise. 

Basic considerations .— Consider the configuration shown in figure k 
where the vertical coordinate of the wing surface is denoted by 
z' = Z(x' ,y',t) (positive downward). The effect of the position and 
motion of the wing may be given by the disturbance— velocity distribution 
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to be superposed on the uniform stream in order to represent the condi- 
tion of tangential flow at the wing surface. This velocity distribution 
normal to the surface (positive upward) is, for small disturbances, 

= i + T ff (1) 

where I is the coordinate in the wind— stream direction. With the use 
of the relation 


BZ = £z_ cbcl Bz_ By* 

Bl Bx ! SI By* SI 


Bz A BZ ‘ . 

t — 7 cos A + t- — - sin A 
Sx* By* 


the vertical velocity at any point is 


w(x* ,y' ,t) 


= U + v 
St 


Sz 

Sx* 


cos A + v 


sin A 
Sy* 


(la) 


Let the wing be bending so that a segment dy* (see fig. 4) is 
displaced from its equilibrium position by an incremental distance h 
(positive down) and also let the wing segment be twisting about the 
elastic axis through an incremental angle 9 (positive leading edge up). 
The position of each point of the segment may be defined, for small 
deflections, by 


Z = h + x'0 (2) 

The velocity distribution normal to the surface, equation (la), conse- 
quently becomes 


w = h + x’0 + v0 cos A + v(a + x*T)sin A 


(3) 


pvy_ 

where a = t— t is the local bending slope of the elastic axis, and is 

oy ^ 

thus analogous to dihedral, and where t = is the local change of 

dy* 

twist of the elastic axis. 


In accordance with assumption (c) the none irculatory— flow velocity 
potentials associated with the vertical— velocity distribution are first 
needed. In equation (3) the terms involving h, 0, and o are constant 
across the chord, whereas those involving 0 and r vary in a linear 
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manner. The none ire ulatory velocity potentials as in reference 6 and 
the new potentials associated with a and t are 


0h = tib\JT- x 2 

9e = v n et\ji - x 2 

0 a = v n CT tan A b^l~— x 2 

00 = ^(f - a)\/r- x 2 

0t = V tan A b 2 (| - a )\fT- x 2 


(M 


where v n = v cos A and x is the nondimens ional chordwise coordinate 
measured from the midchord as in reference 6, related to the coordinate 
x* in the manner 

x = Si + a 
b 

The velocity potential for the circulatory flow associated with the 
wake may be developed on the basis of assumption (c) and the concepts 
for the infinite unswept wing introduced in reference 6. (Thus the 
circulatory— flow pattern for a section dy’ of the finite swept wing 
is to be obtained from the corresponding flow pattern for an infinite 
uniform yawed wing. This infinite wing is assumed to have undergone 
harmonic oscillations for a long time; the full wake is established, 
remains where formed, and consequently is harmonically distributed in 
space. For the infinite uniform yawed wing results for the circulatory 
flow are like those of reference 6 with v replaced by the component 
v n and with the addition of terms to take care of a and t . ) In 
particular, the strength of the wake acting on each section is deter- 
mined by the condition of smooth flow (the velocity remaining finite) 
at the trailing edge. This condition is utilized in the form 

^(0 r + fa) is equal to a finite quantity at the trailing edge; (where 

0 r is the velocity potential due to the vorticity in the wake, and 0 N 
is the total none ire ulatory velocity potential) and leads to a relation 
analogous to equation (YII) of reference 6 involving the basic quantity 

Q = h + v n 0 + v n a tan A + b^- — aj(0* + v n T tan a) which occurs in the 

terms associated with the wake. The net result of these considerations 
is that the circulatory— flow velocity potential may be regarded as 
determined. 
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The pressure difference between upper and lower surfaces of the 
wing at a point x is 



where 0 is in general the total potential, the sum of circulatory— and 
none irculatory— flow potentials. The last term in equation ( 5 ) is the 
product of the component of main— stream velocity taken along the wing 
and the lengthwise change in the velocity potential, and is often neg- 
lected even in steady— flow work. The question of the retention or neg- 
lect of this last term seems partly dependent on the order in which the 
approximations are introduced; specifically, whether velocity potentials 
for the whole flow pattern are found and then the integrated forces are 
determined or whether section forces are first determined and then inte- 
grated. It seems appropriate to retain at least the none irculatory part 
of 0 in the last term of equation (5)- I n view, however, of the 
nature of the approximate treatment of the circulatory potential and of 
the inherent shortcomings of a strip analysis, in particular the neglect 
of lengthwise variations in wake vortex strength, complicating the 
results by also including 0p In this term does not appear worth while. 
(This neglect of 0p and retention of 0n is realized to involve some 
inconsistencies in that account may not be taken of other higher order 
terms associated with lengthwise variation of the wing wake, which may 
be of the same order as terms retained.) Thus equation ( 5 ) becomes 


ebe’ 

For harmonic motion in each degree of freedom, relations for the 
pressure may be integrated over the chord to yield expressions for the 
air forces and moments. For the sake of separating and identifying the 
terms in force and moment expressions which are due solely to the inclu- 
sion of the last term in equation (5 a ) a special bracket is 

employed. Thus these terms may be readily omitted. Numerical checks 
among the calculations made for the present paper showed the effect of 
inclusion of the last term in ( 5 a) on the calculated results to be quite 
small, even for 60 ° of sweepback within the range of other parameters 
investigated. 



cos A + v 


^ 3ln A J 


5y 


(5 a ) 
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The expressions for the aerodynamic lift (positive down) and for 
the moment about the elastic axis (positive leading edge up), each per 
unit length of the wing, are: 

P = — 2rtpv n bcfit + v n 0 + v n a tan A + b^i- — a j(^0 + v n T tan Aj 
itpb2pi + v n 0 + v n cr tan A + jv n a tan A + v^t tan A + 

+ itpb^a 0 + v n f tan A + -[v n f tan A + v n ~ tan 2 /T}-j 

( 6 ) 


v n 2 tan2A T 


M 


a 


= 2irpv n b 2 ^i + aj c|k + v n 0 + v n <j tan A + bQj- - aj[e + v n T tan A^J - 


Jtpv n : 


b3 


[( I - a ) 


0 + — v n r tan A 


+ rtpb3a 


L [ A + v n 


cr tan A + 


[v n 6r tan A + v n 2 T tan A + v n 2 tan 2 yv}J - jtpb^i + a 2 j|l)' 


v n T tan A + -[v n T tan A + v Q 2 |2_ tan2Af- 


( 7 ) 


where 


C = 0(1^) = F(h n ) + iG(l^) 


is the function associated with the wake developed by Theodorsen in 
reference 6; the reduced frequency parameter is defined by 



mb 

v cos A 


(8) 


As has already been stated, the foregoing expressions were developed 
and apply for steady sinusoidal oscillations. 


> 

e = 0 1 (y)e iu)t 


'J 
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The amplitude, velocity, and acceleration in each degree of freedom are 
related as in the h degree of freedom; that ia, 

h = ioii 
h = — <o^h 

Expressions for force and moment.— With the use of such relations 
equations (6) and (7) may he put into the form 

P = -nph 3 cD 2| (hB ch + ©B c0 ) (10) 


Ma » -» pbV ( l > B ah + SB aS ) 


where 


B C h = l A ch + f tan A (- 1 ^)({ + A ch) + 'h V 


B c0 = A ca + 7T b tan A ( A cr) +- 


— tan^A ( — — ^ 

e Sy« ^ 


Bah = ^ Aah + J tan A^-i + A ah^ +- 


b tan 2 A/- a 


h dy* 


V 


B a0 = A aa + W b tan + 




tan 2 a! 




in which the four following coefficients: 


, 1 2G , 2F 

^ch " 1 kn + 1 k . 


■n 


A ca = a + 


A ah = a + 




- 

2 \2 


+ a A, 


l ch 


(11) 
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aa 


1 2 

-8" a 




*n 


i 

f— - a^— _ ( 

1 _ a 2 W - 

1 + a) 2S_ 


\2 K \ 

* K 

V 2 ; k n 2 _ 


are identical with those used in the case of the unswept wing. 
Additionally ^ 


CT 




J 2 - " + w + (l ~ a K3 

- 1 t["(t + {5 + a2 }) + "* a H 


As was previously stated, the special bracket 11 is used to identify 
terms originating in the last term of equation (5a). 

It is of interest to note that equations (6) and (7) reduce, for 
the case of the wing in steady flow (k n = 0), to 


P = — 2itpbv n ‘ 


+ a 


tan A + Tb tan A i- 


+ 2 I 


Jb 8a 

12 


— ^2 ^ — tan 2 a 


tan^ A b 

2 8y’ 


3 


(10a) 


= 2itpb 2 v n 2 |^0 + a tan A)^i + a^ + Tab tan A^jlj a^ 


Jah bo 

tan 2 A — — b 2 (i + a2 

) 1— tan 2 All 

h 

2 \8 / 

1 V JJ 


(11a) 


per unit length of wing. 
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Introduction of modes.— Equations (10) and (11) give the total 
aerodynamic force and moment on a segment of a sweptback wing oscillating 
in a simple harmonic manner. Relations for mechanical equilibrium 
applicable to a wing segment may be set up, but it is preferable to bring 
in directly the three-dimensional-mode considerations. (See for example, 
reference 7.) This end may be readily accomplished by the combined use 
of Rayleigh type approximations and the classical methods of Lagrange. 

The vibrations at flutter are assumed to consist of a combination of 
fixed mode shapes, each mode shape representing a degree of freedom 
associated with a generalized coordinate. The total mechanical energy, 
the potential energy, and the work done by applied forces, aerodynamic 
and structural, are then obtained by the integration of the section 
characteristics over the span. The Rayleigh type approximation enters 
in the representation of the potential energy in terms of the uncoupled 
frequencies . 


As is customary, the modes are introduced into the problem as 
varying sinusoidally with time. For the purpose of simplicity of 
analysis, one bending degree of freedom and one torsion degree of 
freedom are carried through in the present development. Actually, 
any number of degrees of freedom may be added if desired, exactly as 
with an unswept wing. Let the mode shapes be represented by 


h = [fhCy'ilh 
9 = 


( 12 ) 


where h = h 0 e iait is the generalized coordinate in the bending degree 

of freedom, and 0 = 0 o e ia5t is the generalized coordinate in the 

torsion degree of freedom, (in a more general treatment the mode 
shapes must be solved for, but in this procedure fh(y') and. f 0 (y') 
are chosen, ordinarily as real functions of y' . Complex functions 
may be used to represent twisted modes.) The constants h Q and 0 O 
are in general complex and thus signify the phase difference between 
the two degrees of freedom. 

For each degree of freedom an equation of equilibrium may be 
obtained from Lagrange's equation 



( 13 ) 
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where is a generalized coordinate and Q,^ is the corresponding 

generalized force. The kinetic energy of the mechanical system is 


T 



m [ f h(y' )] 2 cLy‘ 


4 - 


ie 2 

2 - 





+ 


• « 

he 



tfh(y’ )] [ f e(y' )] <*7* 


(U) 


where 

m mass of wing per unit length, slugs per foot 

I a mass moment of inertia of wing about its elastic axis 

per unit length, slug— feet^ per foot 

distance of sectional center of gravity from the elastic 
axis, positive rearward, feet 

The potential energy of the mechanical system may be expressed 
in a form not involving bending— torsion cross— stiffness terms: 




(15) 


where 

El bending stiffness , pound— feet^ 

G-J torsional stiffness, pound— feet ^ 

If Rayleigh type approximations are used to introduce frequency, 
the expression for the potential energy may be written in a more 
convenient form: 



mf h 2dy' 



Ve 2d r 


(15a) 
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Another expression for the potential energy is 



(15b) 


The effective spring constants and C a correspond to unit length 

of wing and thus conform to their use in references 5 to 7. The 
constants are effectively defined by 


These effective spring constants are related to the frequencies 
associated with the chosen modes. For so-called uncoupled modes 
the frequencies appropriate to pure modes (obtained by proper 
constraints) are often used. On the other hand, employment of 
the normal or natural modes and frequencies appropriate to them, 
which might be obtained by proper ground test or by calculation, 
may .be preferred. In either case the convenience of not having 
cross— stiffness terms in the potential— energy expression is noted. 

Application is now made to obtain the equation of equilibrium in 
the bending degree of freedom. Equation (13) becomes 


The term represents all the bending forces not derivable from the 

potential— energy function and consists of the aerodynamic forces 
together with the structural damping forces. The virtual work 5W 
done on the wing by these forces as the wing moves through the virtual 
displacements 6h and 50 is: 





(16) 
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6W = 


dy' 


[( P '°hf ») 5h + (% - C a ^ «)“ 

= f ( p ~ ““di 2 — f hi) f h 4 7' & ]i + 

£ (Ma - W* ^ f 0i) f 0 d y’ 5 ® = Q h Sh + % 6£ (17) 

structural damping coefficient for tending vibration 

structural damping coefficient for torsional vibration 

In this expression the aerodynamic forces appropriate to sinusoidal 
oscillations are used. The application of the structural damping as 
in equation ( 17 ) (proportional to deflection and in phase with velocity) 
corresponds to the manner in which it is introduced in reference 5* 

For the half-wing 

<z h = £ ( p - “» 2 ^ f ii) f h -V’ 

£ (£) 3 jj Wh 2 + tf 1 t)( W + Ach ) (tan A)f “ V + 


where 

Sh 

Sx 


= — it c 


_ h — b(tan 2 A)fi, — % 


~ v 2 
ui 


dy' ‘ 


-+ 0 A ca f 0 f h + £A cT b(tan A)f h 0- + 


- 0 


(”i% 5 )' jS(t “ 2A)fh 02 "" 4 « ehfhZ 


dy* 


( 18 ) 


where b r is the semichord at some reference section. Performance of 
the operations indicated in equation ( 16 ) and collection of terms lead 
to the equation of equilibrium in the bending degree of freedom 
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f h 2dy' 
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£e 

dy' 


dy' + 


(19) 


where 

1 _ m 

K v2 

jrph^ 

By a parallel development the equation of equilibrium for the 
torsional degree of freedom may also be obtained 


h 



fa> + 



dy* ‘ 


dy' 




'' -jL.\ 6 tan 2 A A + a 2 


V8 




dy' 2 


► 1 itp'bjAaj 2 = 0 


( 20 ) 
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where r a = y I^mb 2 , (radius of gyration of wing about the elastic 
axis) . 

Determinantal equation for flutter .— Equations ( 19 ) and (20) may be 
rewritten with the use of the nondimens ional coordinate tj = yq-. They 

then are in the form 


(M 2 + 0B 2 )n:pb r 3a32 = 0 

( 19 a) 

(hl>2 + 0.E 2 )npb r ^(n 2 = 0 

(20a) 


where 


a 2 - 


a.o, 2 


- (l) + lg h) J o (£) 

* f °© ! “ *(• t)<« * 
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Bo = V 


— Aca) [jh( ’l )] |j^0 ( Tl ^ 

r lo ' (£)**“ A(A ° t) ^ h<r|) ] ^ dn + 

b 2 pi.O,, .5 . _ . r -1 d 2 F fl 

F-j[ fc) tan A (^ (0] ^ 4n 
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in which F h (t|) = f h (l*Tj) and F 0 (t|) = f 0 (l'T)). 

The borderline condition of flutter, separating damped and 
undamped oscillations, is determined from the nontrivial solution of 
the simultaneous homogeneous equations (19a) and (20a). Such a 
solution corresponds to the fact that mechanical equilibrium exists 
for sinusoidal oscillations at a certain airspeed and with a certain 
frequency. The flutter condition thus is given by the vanishing of 
the determinant of the coefficients 



Application to the case of uniform, cantilever, swept wings is 
made in the next section. 
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Application to Uniform Cantilever Swept Wings 


The first step in the application of the theory is to assume or 
develop the deflection functions to he used. For the purpose of applying 
the analysis to the wing models employed in the experiments it appeared 
reasonable to use for the deflection functions, Fh(q) and Fe(q), 
the uncoupled first bending and first torsion mode shapes of an ideal 
uniform cantilever beam. Although approximations for these mode shapes 
could be used, the analysis utilized the exact expressions developed 
from equations (120) and (I 06 d), respectively, of reference 8 by appli- 
cation of appropriate boundary conditions. 

The bending-mode shape can be written 


where pj = 0 . 5969 * for first bending. The torsion mode shape can be 


where Pg = for first torsion and Cj and Cg are constants. 

The integrals appearing in the determinant elements Ag, Bg, Dg, 
and Eg are 


■ C 1 cosh + - c o s - e l 


sinh Pj + sin P-j_ 


(cos p^t) — cosh P^q) + sinh p^ — sin P^q 


written 


Fg(q) - Cg sin Pgq 


^ 1.0 


F h 2dq = 1. 8554C x 2 


JO 
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6Fq 

drj 
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O.3I83C2 2 



dt) = 1.5926C-, 2 
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The flutter determinant becomes 
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or more conveniently, when columns and rows of the determinant are divided by appropriate 
terms 
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It is interesting to note that the parameters A and Z'/b r appear only 

in the combination in the immediately preceding determinant. 

Z /b r 

The solution of the determinant results in the flutter condition. 


RESULTS AND DISCUSSION 
Experimental Investigation 


Presentation of experimental data. — Results of the experimental 
investigation are listed in detail in tables I to VII, and some 
significant experimental trends are illustrated in figures 5 to 10. 

As a basis for presenting and comparing the test results, the ratio 
of experimental tunnel stream conditions to the reference flutter 
conditions is employed so that the data indicate more clearly combined 
effects of aspect ratio, sweep, and Mach number. As previously 
mentioned, use of the reference flutter speed Vr serves to reduce 
variations in flutter characteristics which arise from changes in 
other parameters, such as density and section properties, which are 
not pertinent to this investigation. (See appendix B.) 

Some effects on flutter speed.- A typical plot showing the effect 
of compressibility on the flutter speed of wings at various angles of 
sweepback is shown in figure 5. These data are from tests of the 
rectangular-plan-form models (type 30) that were swept back by use of 
the rotating mount, for which arrangement the reference flutter speed 
does not vary with either Mach number or sweep angle. Observe the 
large increase in speed ratio at the high sweep angles. 

The data of reference 1 from tests of a rigid, flexibly mounted 
rectangular model having a rotating base are also plotted in figure 5- 
It can be seen that the data from the cantilever models of the present 
paper which had a similar method of sweep are in conformity with the 
data from the flexibly mounted model. This indicates that, for uniform 
wings having the range of parameters involved in these tests, the 
differences due to mode shape are not very great. 

Figure 6 is a cross plot of the data from figure 5 plotted 
against A at a Mach number approximately equal to O.65. The data 
of the swept wings of constant length-to-chord ratio and of the 
sheared swept wings are also included for comparison. The velocity 
ratio V e /Vp is relatively constant at small sweep angles but rises 

noticeably at the large sweep angles. It is pointed out that the 
reference flutter speed Vr may be considered to correspond to a 
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horizontal line at ^ = 1 for the rotated and constant— length— 

V R 

chord— ratio wings, hut for the sheared wings this reference spe ed 
corresponds to a curve decreasing somevdiat less rapidly than ^cos"A as 
A increases. (See appendix B. ) 


The order of magnitude of some three-dimensional effects may 
he noted, from the fact that the shorter wings ^ = 6.5, fig. 6, 

series have higher velocity ratios than the longer wings ^ ^ = 8 . 5 , 

series IV^ . This increase may he due partly to differences in flutter 
modes as well as aerodynamic effects. 


Some effects on flutter frequency .— Figure 7 is a representative 
plot of the flutter-frequency data given in table II. The figure shows 
the variation in flutter— frequency ratio with Mach number for different 
values of sweep angle for the models rotated hack on the special mount. 
The ordinate is the ratio of the experimental flutter frequency to the 
reference flutter frequency f e /fR. It appears that there is a reduc- 
tion in flutter frequency with increase in Mach number and also an 
increase in flutter frequency with increase in sweep. The data from 
reference 1 show the same trend with increase in sweep. Considerably 
more scatter may he noted in the frequency data than in the speed data 
(fig. 5) from the same tests. 


The results of the tests for rotated wings with chordwise lami- 
nations (models kOA, B, C, D) are given in table II. At sweep angles 
up to 30° the values of the speed ratio V e /VR for wings of this 
construction were low (in the neighborhood of 0.9) , and the flutter 
frequency ratios fe/fR were high (of the order of lA) . As these 
results indicate and as visual observation showed, these models 
fluttered in a mode that apparently involved an appreciable propor- 
tion of the second bending mode. The models with spanwise laminations 
(models 30A, B, C, D) also showed indications of this higher flutter 
mode at low sweep angles; however, these models were able to pass 
through the small speed range of higher mode flutter without suffi- 
ciently violent oscillations to cause failure. At a still higher 
speed these models with spanwise laminations fluttered in a lower 
mode resembling a coupling of the torsion and first bending modes. 

This lower mode type of flutter characterized the flutter of both 
the sheared— and constant— length— chord— ratio models. 


For those wing models having the sheared type of balsa construction 
(models 22*, 23, 2h, and 25), the results are more difficult to compare 
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with those of the other models. This difficulty arises chiefly because 
the lightness of the wood produced relatively high mass-density ratios k 
and partly because of the nonhomogeneity of the mixed wood construction. 
For high values of k the flutter— speed coefficient changes rather 
abruptly even for the unswept models (reference 5). The data are never- 
theless included in table I. 

Effect of shift in center-of— gravity position on the flutter speed 
of swept wings .— Results of the investigation of the effects of center- 
of— gravity shift on the flutter speed of swept wings are illustrated in 
figure 8. This figure is a cross plot of the experimental indicated air 
speeds as a function of sweep angle for various center-of-ngravity posi— 


which serves to reduce the scatter resulting from flutter tests at 
different densities of testing medium. The data were taken in the Mach 
number range between 0.14 and 0.44, so that compressibility effects are 
presumably negligible. As in the case of unswept wings, forward movement 
of the center of gravity increases the flutter speed. Again, the flutter 
speed increases with increase in the angle of sweep. 

The models tested at zero sweep angle (models 91 — 91-2* 91 — 3) were 
of different construction from and of larger size than the models tested 
at the higher sweep angles. Because of the manner of plotting the 
results, namely as experimental indicated airspeed (fig. 8) , a compar- 
ison of the results of tests at A = 0° with the results of the tests 
of swept models is not particularly significant. The points at zero 
sweep angle are included, however, to show that the increase in flutter 
speed due to a shift in the center-of-gravity position for the swept 
models is of the same order of magnitude as for the unswept models. For 
the unswept models, the divergence speed and the reference flutter 

speed Vp are fairly near each other, and although the models appeared 

to flutter, the proximity of the flutter speed to the divergence speed 
may have influenced the value of the critical speed. 

The method used to vary the center of gravity (see fig. 1(g)) 
produced two bumps on the airfoil surface. At the low Mach numbers of 
these tests, however, the effect of this roughness on the flutter speed 
is considered negligible. For proper interpretation of figure 8 the 
fact must be kept in mind that the method of varying the location of 
the center of gravity changed the radius of gyration r a and the 

torsional frequency f a . 

T he effect of sweepforward on the critical speed .— An attempt was 
made to determine the variation in flutter speed with angle of sweep- 
forward by testing wings on the mount that could be rotated both back- 
ward and forward. As expected, however, the model tended to diverge at 


tions. The ordinate is the experimental indicated air speed 
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forward sweep angles in spite of the relatively forward position of the 
elastic axis in this D-spar wing. 

Figure 9 shows a plot of the ratio of critical speed to the refer- 
ence flutter speed Vr against sweep angle A. Note the different 
curves for the sweptback and for the sweptforward conditions and the 
sharp reduction in critical speed as the angle of sweepforward is 
increased. The different curves result from two different phenomena. 
When the wing was swept back it fluttered, whereas at forward sweep 
angles it diverged before the flutter speed was reached. Superposed 
on this plot for the negative values of sweep are the results of 
calculations based on an analytical study of divergence (reference 9) . 
Reasonable agreement exists between theory and experiment at forward 
sweep angles. The small difference between the theoretical and experi- 
mental results may perhaps be due to an inaccuracy in determining either 
the position of the elastic axis of the model or the required slope of 
the lift curve or both. 

The divergence speed Vp for the wing at zero sweep angle, as 
calculated by the simplified theory of reference 5, is also plotted in 
figure 9. This calculation is based on the assumption of a two- 
dimensional unswept wing in an incompressible medium. The values of 
the uncoupled torsion frequency and the density of the testing medium 
at time of flutter or divergence are employed. Reference 9 shows that 
a relatively small amount of sweepback raises the divergence speed 
sharply. For convenience, however, the numerical quantity Vp (based 
on the wing at zero sweep) is listed in table I for all the tests. 

Effect of tip modifications .— Tests to investigate some of the 
over— all effects of tip shape were conducted and some results are 
shown in figure 10. Two sweep angles and two length— to— chord ratios 
were used in the experiments conducted at two Mach numbers. It is 
seen that, of the three tip shapes used, namely, tips perpendicular 
to the air stream, perpendicular to the wing leading edge, and parallel 
to the air stream, the wings with tips parallel to the air stream gave 
the highest flutter speeds. 


Discussion and Comparison of Analytical 
and Experimental Results 

Correlation of analytical and experimental results has been made 
for wings swept back in the two different manners; that is, (l) sheared 
back with a constant value of Ag, and (2) rotated back. The two types 
of sheared wings (series I) and two rotated wings (models 3^B and 30D) 
have been analyzed. 
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Results of some solutions of the flutter determinant for a wing 
(model 3C©) on a rotating base at several angles of sweepback are shown 
in figures 11 and 12. Figure 11 shows the flutter— speed coefficient as 
a function of the bending— to— torsion frequency ratio, and figure 12 
shows the flutter frequency ratio as a function of the bending— to— torsion 
frequency ratio. 

The calculated results (for those wings investigated analytically) 
are included in tables I and II. The ratios of experimental to analy- 
tical flutter speeds and flutter frequencies have been plotted against 
the angle of sweep in figures 13 to 1 6 . If an experimental value 
coincides with the corresponding analytically predicted value, the 
ratio will fall at a value of 1.0 on the figures. Deviations of 
experimental results above or below the analytical results appear on 
the figures as ratios respectively greater than or less than 1.0. 

The flutter— speed ratios plotted in figure 13 for the two rotated 
wings show very good agreement between analysis and experiment over 
the range of sweep angle, 0° to 60°. Such good agreement in both the 
trends and in the numerical quantities is gratifying but probably 
should not be expected in general. In view of the discussion of the 
last term in equation (5a) it may be of interest to mention that 
failure to include the terms arising from the last term of equation (5a) 
in the calculations for model 30 B would decrease the ratio Ve/YA 
corresponding to A = 60° by about 3 percent. The flutter frequency 
ratios of figure l4 obtained from the same two rotated wings are in 
good agreement. 

The flutter— speed ratios plotted in figure 15 for the two types 
of sheared wings do not show such good conformity at the low angles 
of sweep, whereas for sweep angles beyond 45 0 the ratios are consider- 
ably nearer to 1.0. The sheared wings are again observed to have a 
constant value of Ag of 2.0 (aspect ratio for the whole wing would 
be 4.0) . For this small value of aspect ratio the finite— span correction 
is appreciable at zero angle of sweep and, if made, would bring better 
agreement at that point. Analysis of the corrections for finite— span 
effects on swept wings requires further consideration. 

Figures 13 and 15 also afford a comparison of the behavior of 
wings swept back in two manners: ( l) rotated back with constant length— 

to— chord ratio but decreasing aspect ratio (fig. 13 ) , and ( 2 ) sheared 
back with constant aspect ratio and increasing length— to— chord ratio 
(fig. 15). A study of these two figures suggests that the length-to- 

chord ratio rather than the aspect ratio ( S — a — may be the relevant 

\Area / 

parameter in determining corrections for finite swept wings. (Admittedly, 
effects of tip shape and root condition are also involved and have not 
been precisely separated.) 
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Figure 16, which refers to the same sheared wings as figure 15, 
shows the ratios of experimental to predicted flutter frequencies. 

The trend is for the ratio to decrease as the angle of sweep increases. 
Table I shows that the flutter frequency fp obtained with Vp 

and used as a reference in a previous section of the paper is not 
significantly different from the frequency f^ predicted by the present 

analysis . 

A few remarks can be made on estimates of ovej>-all trends of the 
flutter speed of swept wings. As a first consideration the conclusion 
may be made that if a rigid infinite yawed wing were mounted on springs 
which permitted it to move vertically as a unit and to rotate about an 

elastic axis, the flutter speed would be proportional to — =— . A 

cos A 

finite yawed wing mounted on similar springs would be expected to have 

a flutter speed lying above the curve of — - — because of finite- 

cos A 

span effects. For a finite sweptback wing clamped at its root, however, 
the greater degree of coupling between bending and torsion adversely 
affects the flutter speed so as to bring the speed below the curve 

of — - — for an infinite wing. This -statement is illustrated in 
cos A 

figure 17 which refers to a wing (model 30B) on a rotating base. The 
ordinate is the ratio of flutter speed at a given angle of sweep to 
the flutter speed calculated at zero angle of sweep. A theoretical 
curve is shown, together with experimentally determined points. Curves 

of — — — and — — are shown for convenience of comparison. The 

cos A Vcos A 

curve for model 30D (not shown in fig. 17) also followed this trend 
quite closely. The foregoing remarks should prove useful for making 
estimates and discussing trends but are not intended to replace more 
complete calculation. In particular, mention may be made, for example, 
that a far-forward location of section center of gravity would lead to 
an entirely different trend. Moreover, as is apparent from the analysis, 
the bending stiffness can play an increasingly significant role with 
increase in the angle of sweep. 

The experiments and calculations deal in general with wings 
having low ratios of natural first bending to first torsion frequencies. 
At high values of the ratio of bending frequency to torsion frequency, 
the position of the elastic axis becomes relatively more significant. 
Additional calculations to develop the theoretical trends are desirable. 
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CONCLUSIONS 


In a discussion and comparison of the results of an investigation 
of the flutter of a group of swept wings, the manner of sweep is 
significant. This paper deals with two main groups of uniform, swept 
wings: rotated wings and sheared wings. In presenting the data, 

employment of a certain reference flutter speed was found convenient. 

The following conclusions seem to apply: 

1. Comparison with experiment indicates that the analysis presented 
is satisfactory for giving the main effects of sweep, at least for nearly 
uniform cantilever wings of moderate length— to— chord ratios. Additional 
calculations are desirable to investigate various theoretical trends. 

2. The coupling between bending and torsion adversely affects the 
flutter speed. The fact, however, that only a part of the forward 
velocity is aerodynamically effective increases the flutter speed. 

Certain approximate relations can be used to estimate some of the 
trends . 

3* Although a precise separation of the effects of Mach number, 
aspect ratio, tip shape, and centei>-of— gravity position has not been 
accomplished, the order of magnitude of some of these combined effects 
has been experimentally determined. Experimental results indicated 
are 


(a) The location of the section center of gravity is an 
important parameter and produces effects for swept wings similar 
to those for unswept wings over the range (30 percent to 70 per- 
cent chord) of locations tested. 

(b) Appreciable differences in flutter speed have been found 
to be due to tip shape. 

(c) The length— to— chord ratio of swept wings is a more 
relevant finite— span parameter than is the aspect ratio. 

(d) Compressibility effects attributable to Mach number are 
fairly small, at least up to a Mach number of 0.8. 

(e) The sweptforward wings could not be made to flutter but 
diverged before the flutter speed was reached. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., September 9 , 19^-8 
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APPENDIX A 

THE EFFECT OF SWEEP ON THE FREQUENCIES OF A CANTILEVER BEAM 


Early In the investigation it was decided to make an experimental 
vibration study of a simple beam at various sweep angles. The uniform, 
plate— like aluminum-^alloy beam shown in figure 18 was used to make the 
study amenable to analysis. Length— to— chord ratios of 6, 3> and 1.5 
were tested, the length Z being defined as the length along the mid- 
chord. A single 60— inch beam was used throughout the investigation, 
the desired length and sweep angle being obtained by clamping the beam 


in the proper position with a 1^ - by 1^- by l4-inch aluminum-alloy 


crossbar. 


Figures 18 and 19 show the variation in modes and frequencies with 
sweep angle. In most cases, an increase in sweep angle increased the 
natural vibration frequencies. As expected, the effect of sweep was 
more pronounced at the smaller values of length— to— chord ratio. The 
fundamental mode was found by striking the beam and measuring the 
frequency with a self— generating vibration pick— up and paper recorder. 
The second and third modes were excited by light-weight electromagnetic 
shakers clamped to the beam. These shakers were attached as close to 
the root as possible to give a node either predominantly spanwise or 
chordwise. The mode with the spanwise node, designated second rode, 
was primarily torsional vibration, whereas the mode with the chordwise 
node, designated third mode, was primarily a second bending vibration. 

The first two bending frequencies and the lowest torsion frequency, 
determined analytically for a straight uniform unswept beam, are plotted 
in figure 19 . Good agreement exists with the experimental results for 
the length— to— chord ratios of 6 and 3> but for a ratio of 1.5 (length 
equal to 12 inches and chord equal to 8 inches) less favorable agreement 
exists. This discrepancy may be attributed to the fact that the beam 
at the short length— to— chord ratio of 1.5 resembled more a plate than 
a beam and did not meet the theoretical assumptions of a perfectly rigid 
base and of simple— beam stress distributions. The data are valid for 
use in comparing the experimental frequencies of the beam when swept 
with the frequencies at zero sweep, which was the purpose of the test. 
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APPENDIX B 


DISCUSSION OF THE REFERENCE FLUTTER SPEED 


For use in comparing data of svept and unswept wings, a reference 
flutter speed Vp> is convenient. This reference flutter speed is the 

flutter speed determined from the simplified theory of reference 5* 
This theory deals with two-dimensional unswept wings in incompress lble 
flow and depends upon a number of wing parameters. The calculations 
in this report utilize parameters of sections perpendicular to the 
leading edge, first bending frequency, uncoupled torsion frequency, 
density of testing medium at time of flutter, and zero damping. 
Symbolically, 


Variation in reference flutter speed with sweep angle for sheared 
swept wings .— The reference flutter speed is independent of sweep angle 
for a homogeneous rotated wing and for homogeneous wings swept back by 
keeping the length— to— chord ratio constant. For a series of homogeneous 
wings swept back by the method of shearing, however, a definite variation 
in reference flutter speed with sweep angle exists since sweeping a wing 
by shearing causes a reduction in chord perpendicular to the wing leading 
edge and an increase in length along the midchord as the angle of sweep 
is increased. The resulting reduction in the mass— density— ratio parameter 
and first bending frequency tends to raise the reference flutter speed, 
whereas the reduction in semichord tends to lower the reference flutter 
speed as the angle of sweep is increased. The final effect upon the 
reference flutter speed depends on the other properites of the wing. 

The purpose of this section is to show the effect of these changes on 
the magnitude of the reference flutter speed for a series of homogeneous 
sheared wings having properties similar to those of the sheared swept 
models used in this paper. 

Let the subscript o refer to properties of the wing at zero 
sweep angle. The following parameters are then functions of the sweep 
angle: 



b = b Q cos A 


cos A 


38 


NACA TN 2121 


Since m is proportional to b. 


irpb . 

k = = K n cos A 

m ° 


Similarly, since I is proportional to b 


S - ^fVir ■- (S) <C0B A)2 


Also, "because f a is independent of A, 


2a » (%A (cos a ) 2 


f \f / 

a Va /_ 


An estimate of the effect on the flutter speed of these changes in 
semichord and mass parameter with sweep angle may be obtained from the 
approximate formula given in reference 5. 


VR 





0.5 

0.5 + a + Xq, 


= Vr qV /cos A 


This approximate analysis of the effect on the reference flutter speed 
does not depend upon the first bending frequency but assumes fh/f a to 

be small. 

In order to include the effect of changes in bending— torsion 
frequency ratio, a more complete analysis must be carried out. Some 
results of a numerical analysis are presented in figure 20, based on a 
homogeneous wing with the following properties at zero sweep angle: 


Xcg 

= 50 


= 0.333 

x ea 

= ^5 

(*) 



= 10 

r 2 

= 0.25 

0 


fa 

= 100 

£). 

= 0.A 


( 
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In figure 20 the curve showing the decrease in Vp with A is slightly 
above the \fc os A factor indicated by the approximate formula. 

Effect of elastic— axis position on reference flutter speed .— As 
pointed out in the definition of elastic axis, the measured locus of 
elastic centers x ea ' fell behind the "section" elastic axis for the 

swept models with bases parallel to the air stream. In order to get 
an idea of the effect of elastic— axis position on the chosen reference 
flutter speed, computations were made both of V R and a second 

reference flutter speed V R ' similar to V R except that x ea ' was 
used in place of x ea - The maximum difference between these two values 

of reference flutter speed was of the order of 7 percent. This differ- 
ence occurred at a sweep angle of 60° when the "wing" elastic axis was 
farthest behind the "section" elastic axis. Thus, for wings of this 
type, the reference flutter speed is not very sensitive to elastic- 
axis position. The reference flutter frequency f R 1 was found in 

conjunction with V R ’ . The maximum difference between f R and f R 1 

was less than 10 percent. Thus, the convenient use of the reference 
flutter speed and reference frequency is not altered by these elastic- 
axis considerations. 


ho 
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TABLE I.- DATA FOE SHEARED SWEPT MODE IE - SERIES I 


Spruce ving8 


Model 


A 

(deg) 


f bl 

(cps) 


**2 
( CDs) 


ft 

(cps) 


(cps) 


GJ 

( lb— in. 2 ) 


El 

(lb-in. 2 ) 


NACA 

airfoil 

section 


l 

(in.) 


(in.) 


b 

(ft) 


x cg 
(percent 
chord) 


(percent 

chord) 


(percent 

chord) 




^SlUfisN 
Vcu ft ) 


Percent 

Freon 


fe 

(cps) 


11A 

11A’ 

1TB’ 

12 

12 

12 

13 

14 

14 

15 
15 


0 

0 

0 

15 

15 

15 

30 

**5 

1+5 

6o 

60 


'1*5 

26 

29 

1*3 

1+2 

1+2 

33 

22 

21 

12 

12 


272 


196 

139 

136 

69 

67 


108 

59 

6l 

103 

105 

103 

94 

93 

92 

93 
93 


107 

37 

1*3 

103 

105 

102 

93 

92 
91 

93 
93 


15,000 


25,100 


14,400 
l4, 400 
14,1+00 
11, 100 
9,240 
9,240 
4,520 
4,520 


54,700 

54,700 

54,700 

53,500 

33,000 

33,000 

19,100 

19,100 


16-005 

16-005 

16-005 

16-005.2 

16-005.2 

16-005.2 

16-005.3 

16-007.1 

16-007.1 

16-010 

16-010 


0.89 

.39 

.39 

.88 

.38 

.88 

.87 

.35 

.35 

.31 

.3i 


16.0 

16.0 

16.0 

16.6 

16.6 

16.6 

13.2 

22.6 
22.6 
32.0 
32.0 


8.0 

3.0 

9.0 
7.72 
7.72 
7.72 
6. 87 
5.62 
5.62 

4.0 
4.0 


0.333 

• 333 

• 333 

• 321 
.321 
.321 
.294 
.234 
.234 

.167 

.167 


48.4 

48.4 

48.4 

43.5 

1+8.5 

1*8.5 

48.3 

48.8 

49.3 

48.9 

48.3 


45 

26.6 

29.7 

46.3 

46.3 

46.3 

46.0 

46.0 

46.0 

46.0 

46.0 


-0.032 

-.032 

-.032 

-.03 

-.03 

-.03 

-.024 

-.024 

-.024 

-.024 

-.024 


- 0.10 

-.468 

-.406 

-.074 

-.074 

-.074 

-.080 

-.080 

-.080 

-.080 

-.090 


O.232 

.396 

• 371 

.23 

.23 

• 23 
•23. 

.23 

.23 

.23 

.23 


13-3 

17.6 

1+0.5 

5.69 

8.47 

11.2 

7.15 

7.78 

19.8 

9.10 

14.0 


0.00287 

.00217 

.000943 

.00725 

.00486 

.OO367 

.00746 

.00720' 

.00285 

.00757 

.00493 


95 
0 

88 

96 

98 

97 

99 
85 
94 
92 
90 


66 

42 

38 

64 

62 

55 

61 

54 

37 

37 

36 


Model 


fR 

(cps) 


fA 

(cp8) 


(deg) 




V e 

(nroh) 


Vr 

(mph) 


V 

(mph) 


(mph) 


V A 


V D 

(mph) 


Remarks 


11A' 

11B' 

12 

12 

12 

13 

14 

14 

15 
15 


70 

40 

40 

70 

71 
69 
60 
56 
51 
53 
51 


71 

69 

65 

61 

53 

58 

55 


0.62 

1.12 

.97 

.63 

.59 

.54 

.66 

.59 

.41 

.40 

.39 


0.9^ 

1.03 

.91 

.92 

.87 

.90 

1.01 

.97 

• 72 

• 70 
.71 


O.87 

.80 

.91* 

.88 

.70 

.64 

.65 


50 

90 

170 

70 

50 

50 

70 

60 

40 

40 

30 


235 

85.O 

70.5 

375 

320 

307 

33b 

300 

234 

265 

264 


0.92 

.24 

• lb 
.64 

• 71 
•19 
.62 
.56 
.31 

• 51 
.62 


274 
191 

262 

218 

245 

276 

202 

196 

275 
179 
222 


260 

129 

197 

176 

206 

225 

154 

13b 

191 

103 

124 


260 

129 

197 


107 

127 


205 

220 

159 

166 

245 

184 


1.80 

3.58 

4.22 

1.54 

1.70 
1.95 
1.77 
2.11 
2.99 

2.71 
3.35 


1.05 

1.48 

1-33 

1.24 

1.19 

1.23 

1.31 

1.46 

1.44 

1.73 

1.79 


1.20 

1.25 

1.27 

1.18 

1.12 

.97 

1.00 


314 

533 

183 

175 

217 

245 

149 

119 

187 

105 

122 . 


Tunnel excitation frequency = 67 cps. 

Model failed. ^ Slotted 4^. inches from trailing edge. 
Model failed. _ Slots uncovered. 

Tunnel excitation frequency = 6l cps. 
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TABLE I.— DATA FOR SHEARED SWEPT MODELS - SERIES I - Concluded 


Balsa wings 








Model 

A 

(deg) 

A g 

f *l 

( CD 8) 

% 

(cos) 

(cps) 

(cps) 

GJ 

(lb— In. 2) 

El 

(lb-in.?) 

NACA 

airfoil 

section 

M C r 

l 

(in.) 

c 

(in. ) 

b 

(ft) 

( percent 
chord) 

x ea 

(percent 

chord) 

x ea* 

(percent 

chord) 

a + x a 

a 

r 2 
r a 

1 

K 

P 

( 8luge\ 

\CU ft/ 

Percent 

Freon 

f e 

(cps) 

22’ 
22' 
22 ' 
23 
23 
23 

23 

24 
24 
24 
24 
24 
25A 
25B 

15 

15 

15 

30 

30 

30 

30 

45 

45 

45 

45 

45 

60 

60 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

31 

31 

31 

35 

34 

34 

34 

19 

19 

19 

19 

19 

8.6 

3.6 

155 

1^4 

154 

219 
216 

220 
216 
123 
122 
122 
120 
120 

54 

48 

63 

64 
64 
89 
39 
91 
39 

73 
75 
75 

74 
73 
66 
70 

61 

62 
62 
39 
89 
91 
89 

73 
75 
75 

74 
73 
65 
68 

6230 

6230 

6230 

6230 

2810 

2310 

2810 

2810 

2810 

1950 

27,900 

27,900 

27,900 

27,900 

10,800 

10.300 
10,800 
10,800 

10.300 
6,470 
5,500 

16-005.2 

16-005.2 

16-005.2 

16-005.3 

16-005.3 

16-005.8 

16-005.8 

16-007.1 

16-007.1 

16-007.1 

16-007.1 

16-007.1 

16-010 

16-010 

0 .89 
.38 
.38 

.87 

.37 

.87 

.37 

.35 

.85 

.85 

.85 

.85 

.81 

.81 

16.6 

16.6 

16.6 

18.2 

18.2 

13.2 

18.2 
21.8 
21.8 
21.3 
21.8 
21.8 
32.0 
32.0 

7.72 
7.72 
7.72 
6.37 
6.37 
6.37 
6.87 
5.66 
5.66 
5.66 
5 .66 
5.66 
4.0 
4.0 

O.321 

.321 

.321 

.234 

.234 

.234 

.234 

.236 

.236 

.236 

.236 

.236 

.167 

.167 

48.8 

48.8 

48.8 
48.0 
48.0 

48.0 
48.0 

47.0 

47.0 

47.0 
47.0 
47.0 

46.9 
46.9 

42.4 

42.4 

42.4 

48.0 

48.0 

48.0 

48.0 

49.0 

49.0 

49.0 

49.0 

49.0 

40.0 
40.0 

42.4 

42.4 

42.4 

52 

52 

52 

52 

57 

57 

57 

57 

57 

71 

71 

-0.024 

-.024 

-.024 

-.04 

-.04 

-.04 

-.04 

-.06 

-.06 

-.06 

-.06 

-.06 

-.062 

-.062 

-O.I52 

”.152 

-.152 

-.04 

-.04 

-.04 

-.04 

-.02 

-.02 

-.02 

-.02 

-.02 

-.20 

-.20 

0.292 

.292 

.292 

.304 

.304 

.304 

.304 

• 311 
.311 
.311 

• 311 
.311 
.359 
.359 

2.19 
3.82 
18.7 
3.18 
8.54 
9.15 
14.9 
3.64 
8.40 
13.2 
29.4 
30.6 
34.6 
‘ 9.36 

0.00854 

.00488 

.00100 

.00864 

.00321 

.00300 

.00184 

.00784 

.00339 

.00216 

.000970 

.000933 

.00095^ 

.OO353 

98 
93 

92 

99 
91 

89 

90 
85 

93 

91 
74 
89 
88 
91 

50 

51 
45 
60 
62 
60 
53 
51 
49 
45 

34 

29 
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TABLE II.- ROTATED WINGS - SERIES II 


Lengthvlse laminations 








Model 

A 

(deg) 

A g 

f hi 

(cps) 

f h2 

(cps) 

ft 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

Mcr 

l 

(in.) 

c 

(in.) 

b 

(ft) 

( recent 

chord) 

x ea 

(percent 

chord) 

^ea 

( percent 
chord) 

a + Xa 

a 

r 2 
r a 

1 

K 

P 

/'slugs \ 
Vcu ft / 

Percent 

Freon 

fe 

(cps) 

30A 

30B 

30B 

3OB 

3OB 

3OB 

3OB 

3OC 

3OC 

3OC 

30C 

3OC 

3OC 

3OC 

3OC 

3OD 

30D 

3OD 

3OD 

30D 

30D 

0 

0 

30 

30 

45 

1+5 

60 

0 

0 

0 

15 

30 

30 

30 

1+5 

15 

15 

15 

30 

1+5 

60 

6.20 
6.20 
1+ .65 
1+.65 
3.10 
3.10 
1.55 
6.20 
6.20 
6.20 
5.78 

I+.65 

I+.65 

1+.65 

3.10 

5.78 

5.78 

5.78 

I+.65 

3.10 

1.55 

11.9 

12.0 

12.1 

12.0 

12.1 

12.2 

12.0 

12.2 

12.2 

13.3 

12.2 

12.2 

12.2 

12.2 

12.2 

13.2 

13.2 

13.2 
13.5 

13.3 
13.5 

76.0 

72.6 

73.0 

73.0 

73.0 

73.0 
72.5 

69.0 

69.0 

70.0 

69.0 

69.0 

70.0 
70.0 
70.0 
80.2 
80.2 
80.2 
81.7 
81.7 
82.0 

90.1+ 
90 >0 

91.0 

90.0 

91.0 

90.0 
90.0 
86.0 
86.0 
81+.0 
86.0 
86.0 
86.5 

86.5 

86.5 
87.1 
87.1 

87.1 

92.5 

88.2 

90.5 

83.0 

88.0 
88.8 
88.0 
88.8 
88.0 
88.0 
75.8 
75.8 
7I+.2 

75.8 

75.8 

76.2 

76.2 

76.2 

82.4 

82.1+ 

82 .4 

87.4 

83.4 

85.5 

3760 

3760 

3760 

376O 

376O 

3760 

3760 

4000 

4000 

4000 

4000 

4000 

4000 

4000 

4000 

4350 

4350 

4350 

4350 

4350 

4350 

6920 

6920 

6920 

6920 

6920 

6920 

6950 

6950 

6950 

6950 

G950 

6950 

6950 

6950 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

0.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

.81 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

24.8 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

0.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

.167 

46.0 

46.0 

46.0 

46.0 

46.0 

46.0 

46.0 

48.5 

48.5 

48.5 

48.5 

48.5 

48.5 

48.5 

48.5 

48 

48 

48 

48 

48 

48 

35 

40 

40 

40 

40 

40 

40 

39 

39 

39 

39 

39 

39 

39 

39 

39.5 

39.5 

39.5 

39.5 

39.5 

39.5 

35 

40 

40 

40 

40 

40 

40 

39 

39 

39 

39 

39 

39 

39 

39 

39.5 

39.5 

39.5 

39.5 

39.5 

39.5 

-0.08 

-.08 

-.08 

-.08 

-.08 

-.08 

-.08 

-.03 

-.03 

-.03 

-.03 

-.03 

-.03 

-.03 

-.03 

-.04 

-.04 

-.04 

-.04 

-.04 

-.04 

— 0 . 30 
-.20 
-.20 
-.20 
-.20 
-.20 
-.20 
-.22 
-.22 
-.22 
-.22 
-.22 
-.22 
-.22 
-.22 
-.21 
-.21 
-.21 
-.21 
-.21 
-.21 

0.311 
.277 
.277 
.277 
.277 
.277 
.277 
.292 
.292 
.292 
. .292 
.292 
.292 
.292 
.292 
.280 
.280 
.280 
.280 
.280 
.280 

36.8 

37.8 

37.7 

37.8 
37.8 

37.8 

39.8 

40.5 

98.9 

92.6 
92.6 

40.0 
81.4 

80.0 
45.2 

8.70 

8.72 

8.76 

8.90 

8.85 

9 . 5 ‘> 

0.00220 

.00214 

.00215 

.00214 

.00214 

.00214 

.00204 

.00200 

.000820 

.000876 

.000870 

.00202 

.OOC995 

.00100 

.00179 

.OO933 

.OO93O 

.00927 

.00910 

.00905 

.00852 

0 

0 

0 

0 

0 

0 

0 

89 

86 

83 

81 

89 

86 

85 

87 

99 

99 

99 

99 

99 

99 

42 

48 

51 

50 

55 

3*+ 

24 

21 

27 

37 

31 

40 

50 

51 
51 
53 
5 6 
65 


Model 

fR 

(cps) 

fA 

(cps) 

fe 

fa 

fe 

fR 

fe 

*A 

<P 

(deg) 

( * b 

M 

V e 

(mph) 

VR 

(mph) 

VR* 

(mph) 

VA 

(mph) 

v e 

hooa 

la 

VR 

la 

TD 

(mph) 

Remarks 

\sq ft; 

3OA 

45 


0.51 

0.91 



70 

127 

0.30 

232 

209 

209 


3.91 

1.11 


318 

Wing failed. Tunnel excitation frequency = 40.7 cps. 

30B 

44 

46 

.54 

1.08 

1.04 

60 

121 

.29 

229 

212 

. 212 

215 

3.64 

1.08 

1.06 

263 


30B 

47 

47 

.57 

1.08 

1.08 

60 

126 

.30 

235 

214 

214 

229 

3.74 

1.10 

1.03 

266 


30B 

44 

47 

.57 

l.l4 

1.08 

40 

120 

.30 

237 

212 

212 

229 

3.77 

1.12 

l.ojf 

263 


3OB 

44 

47 




— 

166 

.34 

269 

214 

214 

265 

4.28 

1.26 

1.01 

266 


30B 

44 

47 

.62 

1.25 

1.16 

50 

169 

i .35 

272 

212 

212 

265 

4.32 

1.28 

1.02 

263 


30B 

46 

48 





275 

! 

1 350 

219 

219 

353 

5.59 

1.60 

•?9 

265 

Wing failed. 

30C 

4l 


.45 

.83 



30 

104 

.63 ! 

' 219 

189 

I89 


4.05 

1.16 

— 

249 


30C 

37 


.32 

.66 

— 

30 

74.4 

.81 

! 286 

290 

290 


5.29 

.986 

— 

393 


30C 

3 6 


.29 

.59 

— 

30 

79.6 

■ .82 , 

288 

270 

270 


5.43 

1.07 

— 

369 

Wing failed. 

30C 

36 


.36 

.74 

— 

30 

72.5 

■ .78 1 

278 

282 

282 


5.13 

.986 

— 

376 


30C 

4l 


.48 

.89 

— 

50 

113 

.65 

226 

187 

187 


4.18 

1.21 

— 

248 


70C 

4l 






— 


88.1 

.81 

284 

263 

263 


5.22 

1.08 

— 

355 


30C 

38 


.40 

.80 

— 

30 

88.6 

.81 

289 

260 

260 


5.32 

i.n 

— 

352 


30C 

4l 


.53 

.98 

— 

40 

147 

.76 

273 

199 

199 


5.02 

1.37 

— 

265 


30D 

51 

51 

.61 

.98 

.98 

50 

no 

.31 

104 

100 

100 

100 

1.77 

1.05 

1.04 

119 


30D 

52 

52 

.61 

.98 

.98 

50 

115 

.32 

107 

100 

100 

101 

1.82 

1.08 

1.06 

119 


30D 

52 

52 

.61 

.98 

.98 

50 

121 

.33 

109 

100 

100 

101 

1.85 

1.10 

1.08 

119 


30D 

54 

56 

.61 

.98 

.95 

40 

150 

.38 

123 

106 

106 

116 

1.97 

1.16 

1.06 

129 


i 30D 

52 

55 

.67 

1.08 

1.02 

60 

178 

.41 

135 

101 

101 

130 

2.26 

1.34 

1.04 

122 


1 30D 

53 

57 

.77 

1.24 

1.14 

90 

307 

.55 

182 

107 

107 

182 

2.98 

1.70 

1.00 

130 



NACA 


uo 


MCA TN 2121 




TABLE II.- ROTATED WINGS - SERIES II - Concluded 


Chordvise laminations 




Model 

A 

(deg) 

A g 

fhi 

(CT>8) 

fh2 

(cps) 

ft 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

Mcr 

Z 

(in.) 

c 

(in.) 

1) 

(ft) 

x cg - 
( percent 
chord) 

x ea 

(percent 

chord) 

x ea* 

(percent 

chord) 

a 4 Xq. 

a 

r a 2 

1 

K 

P 

/slugs \ 
\cu ft/ 

Percent 

Freon 

fe 

(cps) 

40A 

0 

6.20 

9.U 

57. U 

90.0 

88. U 

351*0 

5250 

16-010 

0.81 

2U.8 

U 

0.167 

U6 

UO 

UO 

-0.08 

-0.20 

0.277 

36.5 

0.00222 

0 

62 

UOA 

0 

6.20 

9 6 

57.1 

91.0 

88.5 

351*0 

5250 

16-010 

.81 

2U.8 

U 

.167 

U6 

UO 

UO 

-.08 

-.20 

.277 

2 k . 2 

.OO33U 

90 

56 

UOA 

0 

6.20 

9.6 

57.1 

91.0 

88.5 

351*0 

5250 

16-010 

.81 

2U.8 

U 

.167 

U6 

UO 

UO 

-.08 

-.20 

.277 

37.7 

.00215 

89 

6l 

1 UOA 

0 

6.20 

9.6 

57.1 

91.0 

88.5 

351*0 

5250 

16-010 

.81 

2U.8 

U 

.167 

U6 

UO 

UO 

-.08 

-.20 

.277 

75.0 

.00108 

82 

61 

| UOA 

15 

5.78 

9.3 

55.8 

90.6 

88.2 

351*0 

5250 

16-010 

.81 

2U.8 

U 

.167 

U6 

UO 

UO 

-.08 

-.20 

.277 

35.1 

.00231 

0 

61 

: UOA 

30 

4.65 

9.3 

55.8 

90.6 

88.2 

351*0 

5250 

16-010 

.81 

24. 8 

U 

.167 

U6 

UO 

UO 

-.08 

-.20 

.277 

37.5 

.00216 

0 

— 

; UOB 

0 

6.20 

9.5 

55.0. 

90.5 

85.5 

3T10 

5020 

16-010 

.81 

2U.8 

u 

.167 

U9 

UO 

UO 

-.02 

-.20 

.297 

35.5 

.00228 

0 

61 

j UOC 

0 

6.20 

9.0 

5U.U 

61.0 

5 8.2 

2280 

U350 

16-010 

.81 

2U.8 

u 

.167 

U6 

38.5 

38.5 

-.08 

-.23 

.287 

8.7U 

.00928 

100 

29 

! UOD 

0 

6.20 

9.U 

58.0 

88.9 

8U.0 

3330 

5050 

16-010 

.81 

2U.8 

u 

.167 

U8 

39.5 

39.5 

-,0U 

-.21 

.280 

79.0 

.OOO969 

8U 

62 

| UOD 

15 

5.78 

9.6 

58.3 

88.9 

8U.0 

3330 

5050 

16-010 

.81 

2U.8 

u 

.167 

U8 

39.5 

39.5 

-,0U 

-.21 

.280 

36.2 

.00212 

89 

62 

, UOD 

15 

5.78 

9.5 

57.9 

87.5 

82.6 

3330 

5050 

16-010 

.81 

2U.8 

u 

.167 

U8 

39.5 

39.5 

-.OU 

-.21 

.280 

80.0 

.OOO956 

87 

6l 

1 UOD 

30 

U.65 

9.5 

57.5 

89.0 

8U.1 

3330 

5050 

16-010 

.81 

2U.8 

u 

.167 

U8 

39.5 

39.5 

-.04 

-.21 

.280 

88.2 

.000867 

85 

65 

; uod 

45 

3.10 

9.6 

58.3 

88.9 

8U.0 

333O 

5050 

16-010 

.81 

2U.8 

u 

.167 

U8 

39.5 

39.5 

-,0U 

-.21 

.280 

39.1 

.00196 

86 

32 


Model 

fR 

(cps) 

fA 

(cps) 

£a 

fa 

fe 

fR 

Is 

f A 

<P 

(deg) 

(errs) 

M 

Ve 

(mph) 

Vr 

(mph) 

VR‘ 

(mph) 

V A 

(mph) 

bw a 

la 

VR 

% 

Vd 

(mph) 

Remarks 

UOA 

47 


0.70 

1.33 



lUO 

82.0 

0.24 

188 

211 

211 


2.98 

0.892 


260 


UOA 

49 


.63 

1.15 

-- 

60 

86.7 

.45 

155 

18U 

184 


2.45 

.843 

-- 

212 

Tunnel excitation frequency * 57 cps. 

UOA 

46 


.69 

1.33 

— 

70 

69.2 

.50 

172 

215 

215 


2.72 

.800 

-- 

265 


UOA 

43 



.69 

l.UU 

— 

70 

63.6 

.65 

234 

299 

299 

— 

3.70 

.784 

— 

373 


UOA 

46 


.68 

1.30 

-- 

90 

93.9 

.26 

201 

20 8 

208 


3.19 

.967 

-- 

254 


UOA 

46 

— 

— 

— 

-- 

— 

127 

.30 

• 235 

213 

213 


3.73 

1.10 

-- 

263 

Wing failed. 

UOB 

45 

— 

.71 

1.37 

— 

10 

77.7 

.23 

178 

191 

191 

— 

2.91 

.932 

— 

247 

Wing failed. 

UOC 

36 

— 

.51 

.83 

— 

80 

57.6 

.23 

75.3 

74.5 

74.5 

— 

1.8l 

1.01 

-- 

90.4 

Wing failed. 

UOD 

40 

— 

.73 

1.5U 

— 

30 

52.3 

.62 

221 

281 

281 

— 

3.69 

.787 

-- 

370 

Tunnel excitation frequency = 6l cps. 

UOD 

UU 

— 

.74 

l.Ul 

— 

70 

72.7 

.51 

177 

194 

194 

— 

2.95 

.913 

— 

251 


UOD 

UO 


.74 

1.54 

— 

50 

57-9 

.67 

236 

279 

279 


3.99 

.846 

-- 

367 


UOD 

UO 


.77 

1.63 

-- 

60 

79.4 

.82 

290 

298 

298 


4.83 

.973 

— 

392 


UOD 

UU 


.38 

.73 

— 

80 

138 

.73 

254 

200 

200 


4.24 

1.27 

-- 

261 

Wing failed. 



NACA TN 2121 


TABLE III.- DATA FOR MODELS USED IN SWEEPFORWARD TESTS - SERIES III 


Model 

A 

(deg) 

Ag 

fhi 

(cps) 

f h2 

(cps) 

ft 

( CDS ) 

fa 

(cps) 

GJ 

(lb-in. 2) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

Mcr 

l 

(in.) 

c 

(in.) 

b 

(ft) 

x cg 

(percent 

chord) 

x ea 

(percent 

chord) 

x ea f 
( percent 
chord) 

a + *a 

a 

r 2 
r a 

1 

K 

P 

/ slugs \ 

\CU ft ) 

Percent 

Freon 

fe 

(cps) 

50A 

-30 

4.65 

15 

87 

168 

137 

10,100 

14,100 

16-010 

0.8l 

24.8 

4 

0.167 

50 

33 

33 

0.0 

-0.34 

O.352 

7.98 

0.00895 

96 



50A 

-15 

5.78 

15 

87 

168 

137 

10,100 

14,100 

16-010 

.81 

24.8 

4 

.167 

50 

33 

33 

.0 

-34 

.352 

8.00 

.00892 

96 


50A 

0 

6.20 

15 

87 

163 

133 

10,100 

14,100 

16-010 

.81 

24.8 

4 

.167 

50 

33 

33 

.0 

“34 

.352 

33.1 

.00216 

0 

102 

50B 

0 

6.20 

14 

82 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

.167 

50 

26 

26 

.0 

-.48 

.456 

8.66 

.00823 

99 

91 

50B 

15 

5.78 

14 

80 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

.167 

50 

26 

26 

.0 

-.48 

.456 

8.58 

.00831 

99 

84 

50B 

30 

4.65 

Ik 

80 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

.167 

50 

26 

26 

.0 

-.48 

.456 

9.04 

.00787 

99 

74 

50B 

45 

3.10 

Ik 

80 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

.167 

50 

26 

26 

.0 

-.48 

.456 

9.45 

.00756 

99 

98 


Model 

fR 

(cps) 

fA 

(cps) 

fe 

fa 

fe 

fR 

> 1® 

<P 

(deg) 

Isq ft ] 

M 

v e 

(mph) 

Vr 

(mph) 

Vr’ 

(mph) 

/a , 

(mph) 

v e 

Ve 

Vr 

V e 

v A 

v D 

(mph) 

Remarks 

50A 

98 



.... 

.... 

- 


73.4 

0.26 

86.9 

174 

174 


0.888 

0.498 

- 

294 

Model diverged. 

50A 

98 



.... 

.... 

- 


107 

.31 

105 

174 

174 


1.075 

.603 

- 

294 

Model diverged. 

50A 

79 


0.77 

1.29 

- 

40 

211 

.40 

303 

319 

319 


3.18 

.949 

- 

579 


50B 

94 


.78 

.97 

- 

100 

260 

.52 

170 

172 

172 


2.05 

.989 

- 

704 


50B 

94 


.72 

.90 

- 

70 

257 

.51 

169 

172 

172 


2.04 

.982 

- 

700 

Model failed. 

50B 

93 


.63 

.80 

- 

180 

352 

.61 

202 

179 

179 


2.44 

1.125 

- 

720 


50B 

93 


.84 

1.05 

- 

100 

423 

.68 

226 

179 

179 


2.73 

1.265 

- 

736 




CA TN 2121 


TABLE IV.- SWEPT MODEI£ OF A CONSTANT LENGTE-TO-CHOKD RATIO OF 8.5 - SERIES IV 


£ 


Model 

(deg) 

*8 

f hl 

(cps) 

fh2 

(cps) 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

(lb— in. 2) 

NACA 

airfoil 

section 

Mcr 

Z 

(in.) 

c 

(in.) 

b 

(ft) 

*cg 

(percent 

chord) 

x ea 

( percent 
chord) 

x ea 

(percent 

chord) 

a + Xa 

a 

r 2 

r a 

1 

K 

P 

( slugsN 
\CU ft/ 

Percent 

Freon 

fe 
( cpe 

62 

15 

7.95 

4.9 

29.1 

72.5 

71.8 

3730 

7,820 

16-010 

0.81 

34 

4 

0.167 

4l 

44 

46 

-0.18 

-0.12 

0.175 

13-5 

0.00925 

99 

22 

62 

15 

7.95 

4.9 

29.1 

73.4 

72.5 

3730 

7,820 

16-010 

.81 

34 

4 

.167 

41 

44 

46 

-.18 

-.12 

.175 

37-6 

.00333 

88 

20 

62 

15 

7.95 

4.9 

29.1 

73.4 

72.5 

3730 

7,820 

16-010 

.81 

34 

4 

.167 

4l 

44 

46 

-.18 

-.12 

.175 

59-5 

.00210 

87 

19 

62 

15 

7.95 

4.9 

29.6 

73.5 

72.7 

3730 

7,820 

16-010 

.81 

34 

4 

.167 

4i 

44 

46 

-.18 

-.12 

.175 

130.0 

.000964 

85 

16 

63 

30 

6.38 

4.6 

25.8 

73.5 

73.0 

5450 

5,870 

16-010 

.81 

34 

4 

.167 

4i 

44 

47 

-.18 

-.12 

.175 

15.2 

.00745 

73 

19 

63 

30 

6.38 

3.9 

24.0 

73.0 

72.4 

5450 

5,870 

16-010 

.81 

34 

4 

.167 

4i 

44 

47 

-.18 

-.12 

.175 

26.8 

.00424 

98 

18 

63 

30 

6.38 

4.6 

25.8 

73.5 

73.0 

5450 

5,870 

16-010 

.81 

34 

4 

.167 

4i 

44 

47 

-.18 

-.12 

.175 

46.0 

.00246 

50 

22 

63 

30 

6.38 

4.0 

24.0 

73-0 

72.4 

5450 

5,870 

16-010 

.81 

34 

4 

.167 

4i 

44 

47 

-.18 

-.12 

.175 

53-0 

.00214 

94 

19 

63 

30 

6.38 

4.0 

24.0 

73.0 

72.4 

5450 

5,870 

16-010 

.81 

34 

4 

.167 

4i 

44 

47 

-.18 

-.12 

.175 

98.2 

.00116 

92 

15 

64 

45 

4.75 

4.4 

29.0 

66.0 

65.5 

3500 

6,080 

16-010 

.81 

34 

4 

.167 

4i 

44 

57 

-.18 

-.12 

.175 

50.9 

.00217 

0 

19 

6k 

45 

4.75 

4.2 

27.0 

66.0 

65.5 

3500 

6,080 

16-010 

.81 

34 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

12.1 

.00914 

97 

— 

6k 

45 

4.75 

4.2 

27.0 

66.0 

65.5 

3500 

6,080 

16-010 

.81 

34 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

41.9 

.00263 

54 

18 

6k 

45 

4.75 

4.1 

27.0 

65.0 

64.4 

3500 

6,080 

16-010 

.81 

34 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

51-3 

.00215 

92 

17 

6k 

45 

4.75 

4.1 

27.0 

65.0 

64.4 

3500 

6,080 

16-010 

.81 

34 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

116.0 

.000953 

86 

16 

65 

60 

2.12 

5.7 

33.4 

77.0 

76.2 

4650 

11,980 

16-010 

.81 

34 

4 

.167 

4l 

44 

71 

-.18 

-.12 

.175 

44.1 

.00297 

94 

17 



Model 

fR 

(cps) 

fe 

'fa 

fe 

fR 

<P 

(deg) 

q 

( 1<b ^ 

M 

V e 

(mph) 

Vr 

(mph) 

Vr* 

(mph) 

v e 

bu^ 

Vr 

Vr 

(mph) 

Remarks 

\®q ft; 

62 

35 

0.28 

0.59 

30 

91.8 

0.29 

95.4 

105 

104 

1.85 

0.905 

91.6 


62 

32 

.28 

.64 

20 

73.7 

.41 

143 

167 

171 

2.76 

.856 

153 


62 

31 

.26 

.60 

20 

69.7 

.49 

175 

206 


3.37 

.850 

192 


62 

29 

.22 

.55 

20 

57.5 

.66 

234 

300 


4.50 

.780 

284 


63 

35 

.27 

.56 

180 

98.8 

.29 

111 

111 


2.12 

1.000 

97.6 


63 

33 

.25 

.56 

110 

78.0 

.38 

129 3 

142 


2.49 

.908 

128 


63 

32 

.30 

.69 

180 

82.1 

.40 

176 

183 

— 

3.37 

.962 

170 


63 

31 

.26 

.61 

l40 

74.0 

.52 

179 

195 


3.46 

.918 

180 


63 

29 

.20 

.50 

120 

62.2 

.64 

222 

262 


4.30 

.848 

246 


64 

28 

.29 

.67 

30 

69.6 

.22 

173 

174 

176 

3.69 

.995 

166 


64 

32 




70.6 

.24 

83.9 

91 

90 

1.80 

.923 

81.3 

No record. 

64 

29 

.27 

.61 

0 

68.3 

.36 

155 

160 

160 

3.31 

.988 

132 


64 

27 

.26 

.62 

30 

63.5 

.47 

165 

172 

171 

3.59 

.960 

173 

Record shown in figure 3. 

64 

25 

.25 

.65 

0 

57.5 

.66 

235 

248 


5.10 

.948 

260 


65 

33 

.22 

.51 

0 

172 

.67 

234 

186 

— 

4.29 

1.258 

176 
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TABLE V.- DATA FOR SWEPT MODELS OF A CONSTANT LENGTH-CHORD RATIO OF 6.5 - SERIES V 


Model 

(deg) 

*8 

f hi 

(cps) 

*h2 

(cpe) 

ft 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

Mcr 

Z 

(l-n \ 

c 

t 

(ft) 

x cg 

( percent 
chord) 

x ea 

(percent 

chord) 

x ea f 

( percent 

a + ia 

a 

r a 2 

i 

P 

/slugs \ 

Percent 

Freon 

U 

(cps) 

section 

im,; 

(in.) 

chord) 



K 

t^cu ft J 

72 

15 

6.09 

7.6 

54 

97.3 

96.3 

3730 

7,820 

16-010 

0.8l 

26 

4 

0.167 

4l 

44 

46 

-0.18 

-0.12 

0.175 

37.2 

0.00336 

94 

30 

72 

15 

6.09 

7.6 

54 

97.3 

96.3 

3730 

7,820 

16-010 

.81 

26 

4 

.167 

4l 

44 

46 

-.18 

-.12 

.175 

81.5 

.00153 

89 

22 

73 

30 

4.88 

6.4 

40.0 

98.0 

97.0 

5450 

5,870 

16-010 

.81 

26 

4 

.167 

4l 

44 

47 

-.18 

-.12 

.175 

34.7 

.00327 

96 

29 

73 

30 

4.88 

6.4 

40.0 

98.0 

97.0 

5^50 

5,870 

16-010 

.81 

26 

4 

.167 

4l 

44 

47 

-.18 

-.12 

.175 

57. ^ 

.00198 

95 

24 

73 

30 

4.88 

6.4 

40.0 

98.0 

97.0 

5450 

5,870 

16-010 

.81 

26 

4 

.167 

4l 

44 

47 

-.18 

-.12 

.175 

108 

.00105 

93 

22 

74 

45 

3.25 

6.5 

40.0 

79.0 

78.2 

3500 

6,080 

16-010 

.81 

26 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

14.2 

.00779 

98 

29 

lb 

45 

3.25 

6.7 

39.5 

78.5 

77.7 

3500 

6,080 

16-010 

.81 

26 

4 

.167 

4i 

44 

57 

-.18 

-.12 

.175 

56.0 

.00197 

93 

26 

lb 

b5 

3.25 

6.7 

39.5 

78.5. 

77.7 

3500 

6,080 

16-010 

.81 

26 

4 

.167 

4l 

44 

57 

-.18 

-.12 

.175 

120 

.000923 

90 

21 

75 

60 

1.65 

7.2 

51.8 

82.4 

81.6 

4650 

11,980 

16-010 

.81 

26 

4 

.167 

4i 

44 

71 

-.18 

-.12 

.175 

15.8 

.00829 

95 

39 

75 

60 

1.65 

7.2 

51.8 

84.6 

83.8 

4650 

11,980 

16-010 

.81 

26 

4 

.167 

4l 

44 

71 

-.18 

-.12 

.175 

16.7 

.00783 

100 

39 





Model 

fR 

(cps) 

fe 

fa 

fe 

fR 

'eij 

2. 

(5%) 

M 

v e 

(mph) 

Vr 

(mph) 

Vr* 

(mph) 

7 6 
bu >a 

Je 

^r 

Vd 

(mph) 

Remarks 

72 

43 

0.31 

0.71 

10 

143 

0.59 

197 

220 

221 

2.88 

0.895 

201 


72 

40 

.23 

.55 

0 

109 

.74 

255 

318 

319 

3.73 

.8o4 

297 


73 

^3 

, *3° 

.67 


133 

.57 

193 

216 

214 

2.78 

.893 

196 


73 

4l 

.24 

.57 

80 

118 

.69 

234 

273 


3.38 

.853 

252 


73 

39 

.22 

.55 


90.8 

.82 

280 

363 


4.05 

.770 

3^5 


74 

37 

.37 

.77 

0 

118 

.35 

118 

115 


2.11 

1.025 

111 

Wing failed. 

lb 

33 

.33 

.77 

0 

104 

.64 

219 

214 


3.95 

1.023 

218 


lb 

31 

.28 

.69 

0 

85.5 

.83 

291 

308 


5.24 

.945 

320 


75 

39 

.47 

.99 

30 

294 

.54 

181 

127 

128 

3.11 

1.425 

113 

1 Model damaged at root 

75 

38 

.46 

.97 

0 

295 

.56 

186 

134 

136 

3.05 

1.386 

122 

J Rear half separated from base, 



■f=r 

-J 


NACA TN 2121 


TABLE VI.- DATA FOB TIP-EFFECT MODELS - SERIES VI 


£ 


Model 

A 

(deg) 

A S 

fhl 

(cps) 

fh 2 

(cpe) 

ft 

(cpe) 

fa 

(cpe) 

GJ 

(lb-in. 2 ) 

El 0 
(lb-in. 2 ) 

NACA 

airfoil 

section 

Mcr 

l 

(in.) 

c 

(in.) 

b 

(ft) 

x cg 

( percent 
chord ) 

x ea 

( percent 
chord) 

x ea 

(percent 
chord ) 

a + Xq, 

a 

r a . 2 

1 

K 

P 

/slugsN 
\cu ft/ 

Percent 

Freon 

fe 

(cpb) 

84-1 

45 

3.63 

10 

60 

133 

104 



16-010 

0.81 

29 

4 

0.167 

51 

32 

44 

0.02 

-0.36 

0.378 

9.15 

0.00781 

99 

75 

84-2 

45 

3.63 

10 

61 

135 

107 



16-010 

.81 

29 

4 

.167 

51 

32 

44 

.02 

-.36 

.378 

9.25 

.00764 

99 

60 

84-3 

45 

3.63 

9.6 

58 

118 

93 



16-010 

.81 

29 

4 

.167 

51.5 

32 

44 

.03 

-.36 

.378 

9.55 

.00778 

99 


8^1 

60 

2.75 

5.0 

32 

92 

72 

10,800 

13,400 

16-010 

.81 

44 

4 

.167 

50 

32 

58 

0.0 

-.36 

.378 

34.6 

.00205 

0 

35 

85-2 

60 

2.75 

5.0 

31 

95 

75 

9,850 

12,400 

16-010 

.81 

44 

4 

.167 

' 50 

32 

58 

.0 

-.36 

.378 

34.1 

.00208 

0 

27 

85-3 

60 

2.75 

5.0 

30 

80 

■63 

11,200 

16,600 

16-010 

.81 

44 

4 

.167 

51 

32 

58 

.02 

-.36 

.378 

34.5 

.00207 

0 

22 


Model 

fR 

(cpe) 

fa 

fp 

<P 

(deg) 

(s%) 

M 

V e 

(mph) 

V R 

(mph) 

Vr* 

(mph) 

v e 

tu) a 

^e 

VR 

Vd 

(mph) 

Remarks 

84-1 

76 

0.65 

0.89 

50 

339 

0.60 

199 

142 

— 

2.65 

1.40 

253 

Tip perpendicular to air stream. 
Model failed. 

84-2 

78 

.51 

.70 

0 

382 

.63 

213 

146 


2.80 

1.47 

259 

Tip perpendicular to leading edge. 
Model failed. 

84-3 

68 




346 

.60 

201 

127 


3.02 

1.58 

229 

Tip parallel to air stream. 





Model failed. 

85-1 

43 

.44 

.72 


225 

.41 

322 

185 

189 

6.24 

1.74 

341 

Tip perpendicular to air stream. 
Model failed. 

85-2 

46 

.33 

.54 


173 

.35 

278 

189 

196 

5.21 

1.47 

348 

Tip perpendicular to leading edge. 
Model failed. 

85-3 

?8 

.32 

• 53 

0 

203 

.39 

304 

159 

159 

6.77 

1.91 

295 

Tip parallel to air stream. 
Model failed. 
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TABLE VII.- DATA FOR MODEIS USED TO DETERMINE EFFECT OF CCTTTE3M)F-<SAVITY SHIFT - SERIES VII 


Model 

A 

(deg) 

A g 

?hi 

(cps) 

( cpe) 

ft 

(cps) 

fa 

(cps) 

GJ 

(lb— in. 2 ) 

El 

(lb-in. 2) 

NACA 

airfoil 

section 

l 

(in.) 

c 

(in.) 

(ft) 

( percent 
chord) 

x ea 

( percent 
chord) 

V 1 

x ea 

( percent 
chord) 

a + la, 

a 

*a 2 

1 

K 

P 

/slugs \ 
\cu ft/ 

Percent 

Freon 

fe 

(cps) 

91-1 

91-2 

91-2 

91-2 

91-2 

91-3 

91-3 

91- 3 

92- 1 

0 

0 

0 

0 

0 

0 

0 

0 

15 

6 

6 

6 

6 

6 

6 

6 

6 

6.09 

4.2 
5.5 
5.5 

5.3 
5.5 
5.0 
4.7 
4.7 

8.3 

24 

36 

36 

33 

36 

30 

29 

29 

48 

31 

43 

43 

42 

43 

40 

39 

39 

70 

23 

43 

43 

42 

43 
40 
39 
39 
62 

34,100 

41,200 

41,200 

41,200 

41,200 

28,500 

28,500 

28,500 

3,730 

128,000 
108,300 
108,300 
108,300 
108, 300 
83,700 
83,700 
83,700 

7,820 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

16-010 

Modified 

16-010 

48 

48 

48 

48 

48 

48 

48 

48 

26 

8 

8 

8 

8 

8 

8 

8 

8 

4 

0.333 

.333 

.333 

.333 

.333 

.333 

.333 

.333 

.167 

29.9 

41.0 

41.0 

41.0 

41.0 

49.0 
49.0 
49.0 

31.2 

48 

43.8 

43.8 

43.8 

43.8 

48.4 

48.4 

48.4 

44 

48 

43.8 

43.8 

43.8 

43.8 

48.4 

48.4 

48.4 

46 

-0.402 

-.18 

-.18 

-.18 

-.18 

-.02 

-.02 

-.02 

-.37 6 

-0.04 

-.124 

-124 

-124 

-.124 

-.032 

-.032 

-.032 

-.12 

0.307 

.179 

.179 

.179 

.179 

.160 

.160 

.160 

.298 

17.3 

41.7 

56.4 

12.8 

95.5 

44.3 

36.4 

48.4 

77.9 

0.00871 

.00239 

.00177 

.00783 

.00105 

.00226 

.00274 

.00207 

.00214 

95 

0 

0 

81 

0 

0 

76 

75 

0 

12.5 

16 

16 

20 

15 

18 

15 

14 

26 

92-2 

15 I 

6.09 

8.3 

49 

95 

95 

3,730 

7,820 

Modified 

16-010 

26 

4 

.167 

42.9 

44 

46 

-.142 

-.12 

.136 

76.0 

.00219 

0 

22 

92-3 

15 

6.09 

8.1 

47 

55 

52 

3,730 

7,820 

Modified 

16-010 

26 

4 

.167 

54.5 

44 

46 

.090 

-.12 

.411 

74.5 

.00224 

0 

26 

93-1 

30 

4.42 

6.3 

40 

78 

68 

5,450 

5,870 

Modified 

16-010 

23.6 

4 

.167 

30 

44 

47 

-.40 

-.12 

.310 

78.0 

.00199 

0 

26 

93-2 

30 

4.42 

6.8 

44 

99 

99 

5,450 

5,870 

Modified 

16-010 

23.6 

4 

.167 

43 

44 

47 

-.16 

-.12 

.134 

74.0 

.00210 

0 

23 

93-3 

30 

4.42 

6.3 

51 

54 

50 

5,450 

5,870 

Modified 

16-010 

23.6 

4 

.167 

56 

44 

47 

.12 

-.12 

.428 

73.2 

.00212 

0 

23 

94—1 

-(-45) 

3.81 

4.5 

26 

38 

35 

2,120 

4,520 

Modified 

16-010 

30.5 

4 

.167 

44.5 

56 

— 

-.11 

.12 

.427 

68.2 

.00223 

0 

18 

94-2 

-(-45) 

3.81 

4.8 

28 

70 

70 

2,120 

4,520 

Modified 

16-010 

30.5 

4 

.167 

57.0 

56 

— 

.14 

.12 

.134 

68.2 

.00223 

0 

18 

94-3 

-(-45) 

3.81 

4.6 

28 

40 

38 

2,120 

4,520 

Modified 

16-010 

30.5 

4 

.167 

69.3 

56 

— 

.386 

.12 

.307 

68.2 

.00223 

0 

17 

95'— 1 

60 

1.65 

5.6 

— 

54 

50 

1,900 

4,560 

Modified 

16-010 

26.4 

4 

.167 

31.4 

22 

41 

-.372 

-.56 

.267 

75.8 

.00201 

0 

24 

95’-2 

60 

1.65 

5.9 


71 

47 

1,900 

4,560 

Modified 

16-010 

26.4 

4 

.167 

42.8 

22 

4l 

-.144 

-.56 

.308 

73.0 

.00209 

0 

23 

95*-3 

60 

1.65 

5.8 

35 

40 

27 

1,900 

4,560 

Modified 

16-010 

26.4 

; - 

.167 

54.3 

22 

. 

4l 

.086 

I -.56 

.779 

69.O 

.00218 

0 

23 


I 

Model ! 

( cps) 

ffi 

fa 

£a 

f R 

<P 

(deg) 

( 

\eq ft/' 

M 

^e 

(mph) 

Vr 

(mph) 

Vr' 

(mph) 

v e 

Zft 

% 

(mph) 

Remarks 

91-1 

15 

0.54 

0.82 



153 

0.37 

127 

231 

231 

3.83 

0.548 

79.9 


Model failed. 

91-2 

19 

.37 

.81 

40 

109 

.28 

208 

207 

207 

3.40 

1.000 

192 



91-2 

19 

.38 

.86 

20 

105 

.32 

239 

239 

239 

3.93 

1.000 

224 



91-2 

21 

.47 

.94 

40 

128 

.33 

122 

120 

120 

2.05 

1.02 

104 



91-2 

18 

.35 

.83 

30 

106 

.40 

3O3 

308 

308 

4.97 

.985 

291 



91-3 

17 

.45 

1.09 

100 

61.5 

.20 

159 

158 

158 

2.78 

1.01 

157 



91-3 

17 

.39 

.91 

10 

58.4 

.39 

142 

l4l 

l4l 

2.54 

1.01 

139 



91-3 

16 

• 37 

.89 

0 

57.2 

.44 

163 

161 

161 

2.92 

1.01 

161 



92-1 

36 

.42 

.72 

0 

195 

.38 

293 

415 

422 

6.60 

.706 

245 



92-2 

36 

.23 

.66 

20 

151 

.33 

255 

258 

257 

3.76 

.990 

251 



92-3 

28 

.49 

.93 

20 

87.5 

.25 

191 

176 

177 

5.12 

1.09 

237 



9 3-1 

26 

.39 

.65 


225 

.41 

324 

503 


6.73 

.645 

267 



93-2 

37 

.23 

.64 

70 

156 

.34 

264 

265 


3.72 

.997 

257 



93-3 

27 

.45 

.85 

20 

77.2 

.23 

185 

170 



5.15 

I.09 

231 



94-1 

20 

.51 

.88 

20 

61.0 

.20 

160 

160 



6.38 

1.00 

122 



94-2 

23 

.26 

.78 

— 

62.2 

.21 

162 

139 



3.24 

1.17 

136 


■ Section reversed. 

94-3 

16 

.44 

1.04 

40 

39.5 

.17 

129 

93.2 


4.78 

1.39 

no 



95*— 1 

27 

.49 

.89 

30 

258 

.44 

345 

279 

300 

5.20 

1.24 

00 



95'— 2 

26 

.48 

.86 

20 

212 

.40 

307 

186 

189 

9.15 

1.66 

00 


’ Slotted 2-^ inches from trailing edge. 
16 

95*— 3 

20 

.84 

1.03 

, J°_ . 

125 

.30 

234 

121 

123 

12.1 

1.94 

Li 



NACA 


NACA TN 2121 


VJ1 

o 


Model 15, 25A, 25B 14,24 13, 23 12, 22' 11A, 11A', 118' 



(a) Sheared swept models with a constant geometric aspect ratio of 2. Series I. 
Figure 1.- Model plan form and cross-sectional construction. 
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(b) Models swept back by use of a rotating mount. Series II. 


Figure 1.- Continued. 
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(c) Models in which a rotating mount is used to determine the effect of sweepback and sweepforward 

on the critical velocity. Series III. 


Figure 1.- Continued. 
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(d) Swept models having a length-chord ratio of 8.5. Series IV. 

Figure 1.- Continued. 
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(e) Swept models having a length-chord ratio of 6.5. Series V. 

Figure 1.- Continued. 
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Model 85-1 85-2 85-3 84-1 



(f) Models used to investigate the effect of tip shape on the flutter velocity. Series VI. 
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Figure 1.- Continued. 
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Model A 

(deg) 

91- 1,91-2,91-3* 0 

92- 1,92-2,92-3 15 

93- 1,93.-2,93-3 30 

94- 1,94-2,94-3 45 

95- 1,95-2,95-3 60 

* Chord = 8" lead inside balsa 
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(g) Models used to determine the effect of center -of -gravity shift on the flutter velocity of swept 

wings. Series VII. 


Figure 1.- Concluded. 
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Figure 2.- Model 12 in the tunnel test section, 





Figure 3.- Oscillograph record of model at flutter. 
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Figure 4.- Nonuniform swept wing treated in the present analysis. 
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Figure 5.- Ratio of experimental to reference flutter speed as a function of Mach number for various 
sweep angles for series II models (fig. 1 (b)) on the rotating mount. 
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Figure 6.- Cross plot of ratio of experimental to reference flutter velocity 
as a function of sweep angle for various wings. Mach number is 
approximately 0.65. 





Figure 7 .- Ratio of experimental to reference flutter frequency as a function of Mach number for 
various sweep angles for series II models (fig. 1(b)) on the rotating mount. 
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Figure 8.- Cross plot of flutter speed as a function of sweep angle for several 
center -of -gravity positions. Series VII models (fig. 1(g)). Length-chord 
ratio is approximately 6. 



Comparison of sweepforward and sweepback tests on wings tested on a rotating mount. 

Series III models (fig. 1(c)). 


Figure 9.- 
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Figure 10.- Effect of tip shape on the flutter speed of swept wings. Wings of length-chord ratios 

of 7.25 and 11 (fig. 1(f)). Series VI models. 
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Figure 11.- Theoretical flutter-speed coefficient as a function of 

the ratio of bending to torsion frequency for the rotated model 30B 
at four angles of sweep and with a constant mass -density ratio 


(i = 37-8). 
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Figure 12.- Ratio of theoretical flutter frequency to torsional frequency 
as a function of the ratio of bending to torsion frequency for the 
rotated model 30B at two angles of sweep and with a constant mass- 

density ratio (i = 37 -S)- 
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Figure 13.- Ratio of experimental to theoretically predicted flutter speed as 
a function of sweep angle for two rotated models. 



Figure 14.- Ratio of experimental to theoretically predicted flutter frequency 
as a function of sweep angle for two rotated models. 
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Figure 15.'- Ratio of experimental to theoretically predicted flutter speed as 
a function of sweep angle for two types of sheared models. 



Figure 16.- Ratio of experimental to theoretically predicted flutter frequency 
as a function of sweep angle for two types of sheared wings. 
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Figure 17.- Flutter-speed ratio as a function of sweep angle for model 30 B 

at a constant mass -density ratio = 37.8) , showing analytical and 
experimental results. ^ 
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Node for third mode 
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Figure 18.- Change in nodal lines with sweep and length -chord ratio for the 

• vibration of a dural beam. 
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Figure 19 •- Variation of frequencies with sweep and length-chord ratio 
for the vibration of an aluminum-alloy beam. 
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